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ABSTRACT 


(Gary Huf ford) 

A preliminary design for an unmanned hypersonic research vehicle to 
test scramjet engines is presented. The aircraft will be launched from 
a carrier aircraft at an altitude of 4Q,®9® feet at Mach 9.8. The 
vehicle will then accelerate to Mach 6 at an altitude of 199,99® feet. 
At this stage the prototype scramjet will be employed to accelerate the 
vehicle to Mach 19 and maintain Mach 19 flight for 2 minutes. The 
aircraft will then decelerate and safely land, presumably at NASA Dryden 
Flight Test Center. 


l l 


TABLE OF CONTENTS 


ABSTRACT 

TABLE OF CONTENTS 

LIST OF FIGURES 

INTRODUCTION 

CONFIGURATION 

WEIGHTS ANALYSIS 

PURPOSE 

METHOD 

SYSTEMS 

AERODYNAMIC SURFACES 

BODY STRUCTURE 

THERMAL PROTECTION SYSTEM 

LAUNCH AND LANDING SYSTEM 

PROPULSION SYSTEM 

INLET SYSTEM 

CRYOGENIC FUEL SYSTEMS 

CONTROLS AND AVIONICS 

PROGRAM 

CENTER OF GRAVITY TRAVEL CALCULATIONS 

AERODYNAMICS 

INTRODUCTION 

SUBSONIC 

SUPERSONIC • 

HYPERSONIC 

PROGRAM USE 

LONGITUDINAL STABILITY 

SONIC BOOM 

PROPULSION 

I n t roduc t i on 

Rocke 

Design 

Rocket Design Calculations . . . 

Rocket Spreadsheet Design .... 

Turbo fan/Ramj e t 

Data 

Scramjet . 

Pro file 

Fuel Percentages 


i i 
i i i 


vi 

1 

1 2 

12 

12 

14 

14 

15 

15 

16 
16 
16 
16 
17 
17 
17 

19 

19 

20 
21 
22 

2 3 

29 

32 

34 

34 

3 
3 


39 

40 
42 
46 
46 

i i i 


J <1 CJ> 



GLIDE CALCULATION . 


52 


FLIGHT PLAN 

SCRAM JET INLET DESIGN . 


TURBOFAN /RAM JET INLET . 


60 


COOLING 

Introduction . . 

Cooling Methods 
Conf i gurat i on 


66 

66 

66 

70 


SUPPORT SYSTEMS . 
Control . . 

Commun i ca t inn 


COST ANALYSIS 
CONCLUSION . 
FINAL NOTE 


REFERENCES 


80 


APPENDICES . • . 

Appendix A . 
Appendix B . 
Appendi x C . 1 
Append i x C . 2 
Appendi x C . 3 
Appendix D . 


LIST OF FIGURES 


Figure 1.1 

Figure 1.2 

Figure 1.3; Final Configuration 

Figure 3.1: Current Configuration . . . . 

Figure 3 . 2 : Re f erence Areas 

Figure 5.1: Surface Boon Pressure .... 

Figure 6.2; Rocket Equations 

Figure 6.3: Rocket Calculations 

Figure 6.4: Raw Data (Ttirust/Drag) ... 

Figure 6.5: Fuel Flow Rate 

Figure 6.6: Raw Scramjet Thrust Data . . 

Figure 6.7: Scranjet Fuel Flow Rates . 
Figure 6.8: Turbo fan/Ramje t Spreadsheet . 

Figure 6.9: Altitude vs. Mach No 

Figure 6.10: Rocketed Fuel % 

Figure 6.11: Rocketless Fuel % 

Figure 6.12: Rocketed F vs. M 

Figure 6.13: Rocketless F vs. M 

Figure 6.14: Rocketed Acc vs. M 

Figure 6.15: Rocketless Acc. vs. M . . . 

Figure 6.17: Rocketed M vs. t 

Figure 6.18: Rocketless M vs. t 

Figure 6.19: Alt. vs. Dist 

Figure 6.20: Alt. vs. t 

Figure 6.21: Dist. vs. t 

Figure 6.16: FFR vs. M 

Figure 8.1: Flight Plan 

Figure 10.1: Turbo fan/Ramj e t Inlet . . . 

Figure 10.2: External Ramps 

Figure 10.3: Internal Ramps 

Figure 10.4: Inlet Changes 

Figure 10.5: Turbo f an/Ramj e t Inlet 3-view 

Figure 10.6: Pressure Recovery 

Figure 11.1: Transpiration cooling . . . 

Figure 11.2: Direct cooling 

Figure 11.3: Indirect cooling 

Figure 12.1: Control Diagram 

Figure 13.1: Cost Analysis 

F i gu re A • 1 : Initial configuration . . • « 

Figure A. 2 : Configuration 2 

Figure A. 3: Configuration 3 

Figure A. 4: Configuration 4 

Figure A. 5: Configuration 5 

Figure A. 6: Engine location 

Figure A. 7: Reference values 

Figure A. 8: Cut-away view 

Figure A. 9: Final configuration 


6 

6 

8 

23 

24 
33 

38 

39 
4 1 
4 1 
44 

43 

44 

45 
47 

47 

48 
48 
48 

48 

49 

49 

50 
50 
50 
50 
55 
6 1 
62 
63 

63 

64 

65 
67 

69 

70 
7 5 
76 

83 

84 

85 

86 

87 

88 

89 

90 
9 1 


v 



INTRODUCTION 
(Gary Hufford) 

This report presents a preliminary design for a research vehicle to 
test scramjet engines. The goal of this vehicle is to test scramjet 

engines in a Mach 6 to Mach 10 speed range. 

The purpose of this airplane is to use the vehicle to test scramjet 
engine technology. Since scramjets have not been flight tested, there 
reliability at present is questionable. Many flight tests will be 
necessary in the development of scramjet engines. A Mach 6 test must be 
successful before higher Mach number tests can be conducted since no 
other engine has an operating range above Mach 6. Thus, a vehicle design 
must be constructed that will allow many different Mach number tests 
with a high degree of efficiency. The goal of this design of the 
research vehicle is to minimize the size and weight of the vehicle while 
still allowing the scramjets to be tested in the Mach 6 to Mach 10 speed 
range. The vehicle is to be launched from a carrier aircraft at Mach 
0.8 at 40,000 feet. By minimizing the size of the research vehicle, the 
size of the carrier aircraft can be minimized. A small vehicle (under 
45000 pounds) may even have the capability to be launched from an 
existing airplane such as a B-52. Thus, a small vehicle may greatly 
reduce the cost of a scramjet engine flight test by reducing material 
costs and by possibly eliminating the need to build a new carrier 
aircraft . 

The base design of the research vehicle is given in Figure 1. 
Assuming the engine data given for the turbo f anramj e t. and the scramjet 
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is accurate, this vehicle has a maximum Mach number of 9.35. A two 
minute cruise can be conducted at any Mach number less than or equal to 
Mach 8.75. This design will allow initial scramjet tests to be 

conduc ted . 

Once the scramjets have been proven effective at Mach 8.75, rocket 
attachments can be made to the base design of Figure 1.2. These attach- 
ments are shown in Figure 1.3. The purpose of the rocket attachments is 
to increase the thrust in the transonic flight regime. An increase in 
thrust will eliminate the thrust pinch normally associated with this 
flight regime. A large increase in acceleration will reduce the fuel 
used by the turbof anram je t during transonic flight and will allow more 
fuel to be used by the scramjets at higher Mach numbers. 

The rocket attachments will increase the aerodynamic drag and the 
weight of the base design of the research vehicle. Thus, it is 
advantageous to eject the rockets after their fuel has been exhausted. 
Therefore, the drag and weight increases will be limited to the region 
where the rockets are in operation. Thus, the thrust gain will far 
outweigh the drag and weight penalty during the rockets region of 
operation. After the rocket fuel has been exhausted and the rockets 
ejected, the drag and weight of the vehicle will be the same as the base 
design of the research vehicle. 

The maximum Mach number of the design with rocket attachment is 
9.9. The maximum Mach number for a two minute cruise is 9.25. Thus, the 
upper speed range for a scramjet test is increased by Mach .5 for the 
same vehicle design. Mach 9.9 is the maximum value the scramjet may be 
tested for any period of time for this vehicle. 
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To obtain higher Mach number tests another scramjet engine must be 
added to the design (presently we have three scramjet engines in the 
design). However, our airplane is not large enough to place another 
scramjet engine on the top of the aircraft. Another engine would reduce 
lifting surface area to unacceptable levels. Therefore, the aircraft 
must be enlarged or redesigned to add another engine. This report 
focuses on the base design and base design with rocket attachments. 
Topics included in this report include the determination of the aircraft 
configuration, weight determination, cost analysis, aerodynamic 
analysis, stability analysis, propulsion system analysis, mission 
profile determination, inlet design, cooling considerations, sonic boom 
analysis and internal systems. 
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CONFIGURATION 
(Kevin Draper) 

During our initial design phase, many different configurations were 
considered. Information on delta wing, wing body, and lifting body 
arrangements was studied. Initial drag estimates, in conjunction with 
the General Electric Turbof anramjet engine data, directed our design 
toward incorporating two full scale G.E. Turbof anramj e t engines and one 
SCRAMJET engine. Since the SCRAMJET engine information hadn’t been 
released at this point, it was assumed that the SCRAMJET engine would be 
no larger than the G.E. Turbof anramjet engines being used. Initially, 
a combination deltawing and lifting body configuration with a leading 
edge sweep of seventy-five degrees was the overall design objective. 

As our design configuration was drawn to scale, a major problem was 
encountered -- the plane’s fuselage would have to be very large just to 
hold the full scale engines. The engine configuration being considered 
contained all three engines, side by side, on the underside of the 
aircraft. When drawn to scale, the engines -- without any inlet or 
isentropic ramp -- would barely fit within the delta wing configuration. 
At this point, alternate propulsion systems and arrangements were 
considered. Solid and liquid rockets were examined as alternatives. 
4fter research, it was decided that they were much too heavy to be a 
feasible alternative. Further information concerning a rocket propul- 
sion system is presented later on in this paper. 

Next, it was discovered that the initial drag estimates had been 
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too high. Since the idea of using a single inlet and diverting the 
airflow to the desired engines was still one of our design 
objectives, we decided to keep the twin turbo f anramj e t layout. The main 
reason for keeping this layout was the fact that it permitted the 
engines to use a single air intake. With twin engines, a symmetric 
thrust could be achieved for the turbof anramjet phase of the flight as 
well as in the SCRAMJET flight regime. With these concepts, and the 
lower drag, it was decided to use two G.E. turbof anramje t engines scaled 
to sixty percent size. The scaled engines would fit within the 
configuration that we had been wanting to use. 

When the turbof anramjet engines were scaled, so too was the allowed 
size for the SCRAMJET engine. This allowed the frontal cross section of 
the aircraft to be smaller. These two preliminary designs can be seen 
in figures A.l and A. 2. The planform changes that were made between 
these two versions increased both the internal volume and the planform 
area. The planform area was increased in an effort to improve the 
aerodynamic characteristics of the design. The volume increase was much 
needed due to the voluminous nature of liquid hydrogen fuel. 

Figure A. 3 contains the altered version that resulted once the data 
for the SCRAMJET engine became available. In this version, the 
frontal area of the aircraft was further reduced, lowering drag 
once more. The isentropic ramp was extended beyond the trailing 
edge of the wing to avoid increasing wing planform area. This 
prevented the aspect ratio from becoming sma ller than two. Again, there 
was a serious problem with incorporating three engines, cooling system, 
airframe structure, control system, and the required storage for the 
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hydrogen fuel. As the design process progressed, the aircraft kept 
getting larger, even after the engines were scaled down. 

Instead of continuing the "airplane growth" cycle that had taken 
over our design efforts, it was decided to alter our design and 
investigate new engine arrangements. The over-under engine 
arrangement that had been decided against in the beginning became 
more appealing. Initially, the over-under engine concept had been 
abandoned for two main reasons. First of all, this arrangement would 
require two separate inlets and would require separate mechanisms to 
seal off the various inlets during the different phases of the flight. 
The second reason that this arrangement was not pursued initially was 
the concern for stability and control. With the over-under arrangement, 
it would not have been possible to design the nose of the aircraft to 
assist in balancing the pitching moments from the engines. However, 
once this over-under arrangement was given more consideration, it 
appeared more practical than our previous configuration. The SC'RAMJET 
engine could be mounted on the top of the aircraft and the top portion 
of the fuselage could then be angled to maximize the inlet for the 
SCRAMJET engine’s Mach number range. Similarly, the bottom portion of 
the fuselage could be angled to obtain the desired shock angle for the 
lower mach number ranges under which the turbo f anramj e t engine would be 
in operation. In order to maximize the thrust performance of both 

engines, an isentropic ramp needed to be incorporated into the design. 
The isentropic ramp would contribute to balancing the pitching moments 
since the expanding exhaust would provide a downward force on the tail 
of the aircraft. This downward force would help counteract the nose 
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down pitching caused by having the engine mounted above the airplane’s 
center of gravity. If the aircraft’s pitching behavior needs further 
assistance, the engine could be tilted slightly to lower its moment arm 
with the center of gravity. The separate inlets will cause a signifi- 
cant amount of drag during the flight phases in which they are not in 
use. However, the decision to alter our engine configuration also 
allowed us to use one G.E. turbo fanramje t at ninety percent size instead 
of the twin engine configuration. This allowed a reduction in the 
overall intake area. With these changes, the new design’s frontal area 
was not altered significantly from the old configuration. The overall 
size of the aircraft remained essentially the same, but the internal 
volume of the aircraft was increased. The reduction in the number of 
engines freed up a considerable amount of space within the aircraft. 
Also, the wing was contoured into the fuselage making it more of a 
lifting body than had been previously planned. This allowed an 
elliptical aft fuselage which increased the aircraft's internal volume. 

This new design initially incorporated four angled control surfaces 
which served the multiple role of a lifting surface as well as horizon- 
tal and vertical control surfaces. The lower angled surfaces were 
designed to extend low enough so that skids could be mounted on them for 
the landing of the aircraft. This version of the design is shown at the 
end of this section. The use of these surfaces as skids during landing 
would have required a very strong wingtip structure and would have 
increased the aircraft weight. In addition to this, the control surface 
tips for this arrangement would have extended outside ot the three 
dimensional shock cone. For these reasons, the lower control surfaces 
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were removed and the wing was extended and contoured into the fuselage. 
Also the wing wasn’t allowed to project outside of the shock cone since 
it would not have remained with the aircraft for long if this had been 
allowed to happen! With this design, the skids will have to be mounted 
in the lower surface of the fuselage. The upper control surfaces were 
also shortened to ensure that they would remain inside of the three 
dimensional shock cone. 

Various schematics of our current design are shown in figures 
A. 5, A. 6, A. 7, A. 8, and A. 9. Figure A. 5 gives the overall design 
configuration along with our weight estimate. Figure A. 6 shows the same 
information as figure A. 5 but also shows the internal location of the 
engines in all three views. In figure A. 7, the distinction made between 
the aircraft’s fuselage and wing is shown. This distinction was made 
for the purposes of aerodynamic analysis and much of the information 
calculated based on aircraft geometry is also given in this figure. 
Figure A. 8 contains a cross-section down the center of the aircraft in 
order to further illustrate the engine arrangement as well as the 
location of the forward landing gear. Figure A. 9 contains information 
resulting from our final aerodynamic and propulsive analysis. The 
£» 0 nter of gravitv location is also shown in this figure. 

As the drag, propulsion, and fuel storage calculations were 
refined, it became evident that there still wasn't sufficient fuel 
storage. The initial response to this situation was to modify the 
aircraft fuselage to gain as much volume as possible without making the 
vehicle excessively large. Two larger versions, shown in fi gu res 1.1 
and 1.2. were considered. It became obvious that the design was heading 
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back into the same "growth" cycle that had been encountered with our 
earlier designs. Since one of the goals in the design was to keep the 
aircraft as small and light weight as possible, an alternate approach 
had to be found. 

The design of the propulsion system and the approach used in 
determining the required fuel was altered. Instead of calculating the 
required fuel to fly certain trajectories and maintain test conditions, 
the available fuel volume was used as a constraint. Since fuel volume 
was critical and the SCRAMJET Mach 10 cruise range could not be met, it 
was decided to add solid rockets to the design. This design is also 
flexible and would facilitate preliminary testing of the SCRAMJET 
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engines at lower Mach numbers since the fuel storage is sufficient to 
reach the low end of the testing regime without solid rockets. During 
later stages of development solid rockets can be added to achieve 
testing at higher speed ranges. The details of the propulsion system 
and analysis is given in the section on the propulsion system. The end 
result of this new approach was that the aircraft would be kept at a 
length of 6Q feet and an available fuel volume of l,5QQ cubic feet. 

The current design configuration is given in figures 1.3 and the 
figures at the end of the section. Also shown is the aircraft version 
with the solid rockets, as well as various performance parameters for 

The aircraft version without the solid rockets as 



this configuration. 




Figure 1.3: Final Configuration 


well as the performance parameters for this version, and the fuel 
storage system and preliminary landing skids/wheel locations are given. 
The internal engine locations are also shown. 

The current design still has volume limitations but it is a much 
more versatile design and adaptations can be made by altering the size 
of the solid rockets being used. The volume available with this design 
is fairly close to being a workable configuration. The limiting 
conflict at this stage is being able to accommodate both the control 
equipment as well as the extensive cooling system that is required by 
such a high speed aircraft. 
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Chapter 2 


WEIGHTS ANALYSIS 

(Scott Dallinga) 

PURPOSE: 

In the second design phase, the weight analysis tasks were divided 
into three areas of study. Component weights of eight major systems 
were calculated for each updated conf igurat ion ( shown in table 8.1). 
Also, the position of each component in the aircraft was estimated to 
locate the center of gravity. Additionally, the center of gravity 
travel was calculated with fuel burn to verify the fuel tank configura- 
t i on . 

METHOD : 

Calculation of most component weights were accomplished using a 
program based on WAATS , Weight Analysis for Advanced Transportation 
Systems, by C.R. Glatt. The original code is described in NASA 
Contractor Report CR-242Q. The WAATS method uses statistical analysis 
of data from existing aircraft to calculate component weights of similar 
aircraft configurations. Correlating parameters are based on geometric 
characteristics and aircraft weight. Since the final weight is not 
known, an iterative process is used to calculate the aircraft gross 
weight. Power law relations are used to relate predominant physical 
characteristics affecting the weights of the components. A power law 
has an empirical constant multiplied by the relevant parameter value 
raised to an empirically defined exponent. 



WEIGHT STATEMENT 
Gross Takeoff Weight: 
Empty Weight : 

Gross Fuel Weight: 

Useable Fuel Weight: 

COMPONENT 

AERODYNAMIC 
Wing Structure 
Vertical Stabilizers 

STRUCTURE 

Main Body Structure 
Secondary Structure 
Thrust Structure 

THERMAL PROTECTION 
Structural Insulation 
Body Cover Panels 
Active Cooling System 

LAUNCH AND LANDING GEAR 
Launch Gear 
Landing Gear 

PROPULSION 

1 Turbofanramjet 

3 SCRAM JET Modules 

2 Solid Rocket Boosters 

INLET SYSTEM 
Variable Ramps 
Internal Ducts 

CRYOGENIC FUEL SYSTEM 
Fuel Distribution System 
Tank System A 
Tank System B 
Tank System C' 

CONTROLS AND AVIONICS 
Aerodynamic Controls 
Elect ron i c s 
Hydraulic System 
Av i on i c s 


44 ,732 
30 , 590 
7,143 
7,180 

WEI GHT 


885 

1,560 


4,307 
9 1 5 
56 1 


1,050 

2,150 

2,000 


1 12 
1 , 008 


5,400 

4,740 

7,000 


73 1 
592 


679 
672 
4 10 
500 


572 
1,919 
340 
7 17 



SYSTEMS : 


The original WAATS program was modified to be used with the 
hypersonic research aircraft concept. There are eight major component 
areas used to obtain a weight bui 1 dup * 


1. Aerodynamic surfaces 

2. Body structure 

3. Thermal protection systems 

4. Launch and landing gear 

5 . Engine sys tem 

6. Inlet system 

7. Cryogenic fuel system 

8. Controls and avionics 


AERODYNAMIC SURFACES: 

Part of the problem in calculating the wing weight for the blue 
group’s hypersonic vehicle is that there is no real distinction between 
the body and the wing. To properly apply the power law relations for 
wing weight, the wing area was defined as the outer thin region of the 
body where the thickness is less than seven and one-half inches. The 
bending moment at the wing root is the determining factor for wing 
weight. To properly apply the empirical relations, the wing span for 
structural considerations is two times the width of the body from the 
point where the thickness is 7.5 inches outboard to the edge. The 
coefficients used in the power law relation are for high-speed, 
h i gh- t emperature aircraft such as the X-15, F-1Q8A, and the XB-7Q. 
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BODY STRUCTURE: 


The three components which make up the body structure are the main 
structure, the secondary structure and the thrust structure. The engine 
compartment, avionics bay, and fuel tank area constitute the main body 
area of the plane. Secondary structures include the afterbody, which 
will be used as an expansion surface for the propulsion system. Landing 
gear doors are also included as secondary structures. The secondary 
structure weight is related to the size of these parts. The engines are 
supported by the thrust structure, whose weight is a function of the 
maximum thrust of the engines. 

THERMAL PROTECTION SYSTEM: 

A common way to keep the structural temperature within acceptable 
limits to use insulation with radiative cover panels installed to 
protect the insulation from flight conditions. Most likely, large 
portions of any Mach 10 aircraft will have to be actively cooled, 
perhaps in combination with the use of insulation and cover panels. At 
this point, due to the lack of detailed thermal analysis, an arbitrary 
value for the area covered by the TPS was selected. To be sure, the 
entire underbody of the vehicle will have to be protected from the high 
temperatures caused by air compression. The afterbody will need cooling 
to protect it from the hot gases ejected by the propulsion systems. 
Included with the thermal protection system is an allowance for the 
active cooling system. Two thousand pounds were allowed for this 
sys tern. 
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LAUNCH AND LANDING SYSTEM: 

In our case, landing gear weight is a function of landing weight or 
launch weight minus fuel weight. Launch and landing gear weights are 
related to similar vehicles which are air-launched. 

PROPULSION SYSTEM: 

The projected weights of the turbofanramjets and scramjets have 
been provided by General Electric. Engine mounts are sized according to 
maximum engine thrust. Two solid rocket boosters with a total weight of 
seven thousand pounds are included in the propulsion system weight. The 
boosters count towards the gross takeoff weight but are not included in 
the empty weight since they are ejected as they burn out. 

INLET SYSTEM: 

Average values were selected for the turbofanramjet inlet parame- 
ters since these parameters change as the aircraft accelerates through 
a range of Mach numbers and altitudes. An two dimensional inlet design 
was selected for the turbofanramjet engine. The scramjet inlet weight 
estimate was also based on a 2-D design with parameters similar to the 
turbofanramjet inlet. 

CRYOGENIC FUEL SYSTEMS: 

Included in the cryogenic fuel system are fuel tanks, fuel tank 
insulation and the fuel distribution system. Fuel tank weight is based 
on a three-tank system with a total volume of 1,5Q9 cubic feet. The 
weight is a function of fuel tank volume. The insulation weight is a 
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function of fuel tank wetted area and fuel tank radiating temperature. 
The fuel tanks were shaped to fit in the available area of the plane. 
The fuel tank radiating temperature was estimated at 5QQ degrees F. 
since the structure is thermally protected, the fuel tank should not 
experience a value too much higher than this. An empirical correction 
factor was applied to the insulation weight since a flight duration of 
only 30 minutes is expected. A percentage of the useable fuel weight 
was allowed for liquid hydrogen boiloff. 

CONTROLS AND AVIONICS: 

Values for the weights of avionics, hydraulics, electrical, and 
aerodynamic orientation control systems were calculated based on 
empirical constants used in power law relations for similar aircraft. 
Aircraft size was the most important parameter involved in the empirical 
weight relations for these systems. 

PROGRAM : 

Power- law relations for components that were used in our design 
were taken from the WAATS program. For each weight estimate a data file 
containing relevant geometrical data as well as coefficients and 
exponents for power law relations was read into the derived weight 
analysis program. The data file contains the parameter values as well 
as a description of each variable and the variable’s name used in the 
program to help future program users understand the methodology used. 
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CENTER OF GRAVITY TRAVEL CALCULATIONS: 


To verify the fuel tank configuration of the aircraft, the center 
of gravity was determined with varying amounts of fuel in the tank 
system. Since a forward C.G. is desirable for stability considerations, 
the order of tank drainage was rear to front. The tanks were drained in 
the following order; A, C, then B. The rocket boosters were ejected 
after 200 lbs of fuel were burned. The rocket assisted phase of the 
flight experienced the farthest C.G. shift due to the high weight loss 
and aft position of the seven thousand pound boosters. C.G. shifts due 
to liquid hydrogen fuel burn were found fairly insignificant. Only 
about 19 percent of the total weight of the plane without the solid 
rocket boosters consists of liquid hydrogen. 
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Chapter 3 


AERODYNAMICS 
(John Aslan and Gar/ Huf ford) 


INTRODUCTION 

In order to determine the preliminary aerodynamic performance of 
our vehicle, we developed three separate computer codes, one each for 
the subsonic, supersonic and hypersonic flight regimes. The use of the 
codes allow the user to change any geometrical aspect of the flight 
vehicle and quickly get the performance data for further study. Data 
from these codes can also be compiled and printed in graphical form. 
The codes were developed from various theories and will be described 
further . 

There are a few common points in all three of the codes. First, 
the standard atmosphere tables are generated at 599 foot intervals. 
This is done by finding the temperature of each geometric height in the 
atmosphere. The gradient is either isothermal or it varies linearly 
with height. The altitudes at which changes occur can be found in any 
book where the standard atmosphere tables are presented. There are 
three separate regions of variation between 9 and 199999 feet, two vary 
linearly and one is isothermal. Once the temperatures are found, the 
ratio of a pressure at that specific altitude to a reference pressure at 
a reference altitude is found. This is calculated from an integrated 
form of the hydrostatic equation, depending only on temperature and 
geometric altitude. Then, the thermally perfect gas law is applied to 
calculate the density. These are generated in 3 separate parts in the 


19 


code corresponding to the three different regions of temperature 
gradients . 

Next, the dynamic pressures at every altitude between Q and 1 QOQQQ 
feet in 5QQ foot increments and the Mach numbers between Q and ten in 
0.1 increments are found. The geometry of the vehicle is now read into 
the code. The user is prompted to enter a range of altitudes, Mach 
numbers and angles of attack. The user is then prompted to choose the 
maximum drag and lift and the minimum lift to cut out the data not 
pertinent to the particular aspect of the performance calculations. At 
the end of the program, various data such as lift coefficient, drag 
coefficient etc. are written to separate files for plotting purposes . 
Also, an overview including all data is written to a file for visual 
inspection of the relationship between all the data in that particular 

ins tance . 

SUBSONIC 

The subsonic theories used are compiled in ’'Fundamentals of 
Aircraft Design" by Leland M. Nicolai, 1975. Portions of the theory are 
presented in graphical form only, so the computer prompts the user for 
these data. Also, it should be mentioned that the aircraft was modeled 
as a delta wing configuration with thickness, but no camber. The 

resultant program is in appendix C.l. 

For the lift calculations, linear two-dimensional theory corrected 
for wing sweep, compressibility and finite span were used. Vving-body 
interference effects are taken into account and are dependent on span, 
body diameter and Mach number. The lift interference factor is 
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determined from these and presented in graphical form in the reference. 
The user is prompted for the lift interference factor by the code. 
Since the aircraft design is of low aspect ratio, non-linear lift 
effects are also taken into account. This is done using non-linear 
theory, dependent on planform shape and aspect ratio, presented in 
graphical form in the reference. The computer prompts the user for this 
information. From this data, a lift-curve slope is found. 

The flow over the entire wing was assumed to be turbulent, so the 
Prandtl approximation for the skin friction coefficient due to turbu- 
lence was used. The wing and body zero-lift drag coefficients are 
calculated separately. The Reynolds numbers were based on the mean 
aerodynamic chord and the body length respectively. The remainder of 
the drag is of that due to lift. For this preliminary design, the 
viscous and inviscid effects were combined into one parameter, dependent 
on the wing efficiency factor and aspect ratio. The computer prompts 
the user for the wing efficiency factor, which is presented in graphical 
form in the reference. The sum of the zero-lift drag coefficient and 
the drag coefficient due to lift give the total drag coefficient. 

SUPERSONIC 

The supersonic theories used are also compiled in "Fundamentals of 
Aircraft Design" by Leland M. Nicolai, 1975. Portions of the theory are 
also only presented in graphical form and must be entered by the user of 
the code. As previously indicated, the aircraft was modeled as a delta 
wing configuration with thickness, but no camber. This code is in 
append i x C . 2 . 
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The supersonic lift coefficient is based on supersonic linear 
theory, corrected for three-dimensional flow effects. Theoretical 
normal force coefficients dependent on aspect ratio, sweep and the 
supersonic similarity parameter are presented graphically in the 
reference. The user is prompted for these values by the computer, which 
is done at every Mach number. From this, a supersonic lift-curve slope 
is generated at every Mach number* 

The drag coefficient is also modeled after supersonic linear 
theory. The wing leading edge was chosen to be more sharp than round, 
so the shock wave is assumed attached. Furthermore, the wing is also 
modeled as having a double-parabolic arc section. The computer next 
calculates whether the wing leading edge sees supersonic or subsonic 
flow, which depends on sweep. Once this is determined, the drag 
coefficient due to the shock wave is determined. The skin friction drag 
coefficient for the wing and body are then determined, in the same 
fashion described in the subsonic section. 

Next, the drag coefficient due to lift is found. For this 
preliminary design, the viscous and inviscid effects are combined into 
one parameter, which is dependent on the lift-curve slope. The 
zero-lift drag coefficient and the drag coefficient due to lift are 
added together to give the total drag coefficient. 


HYPERSONIC 

The hypersonic code, listed in appendix C.3, was written from 
theory presented in "Hypersonic Aerodynamics" by R.Truitt, 1959. It is 
extracted from Newtonian flow as a two-dimensional, three-term series 
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foil. Both lift and drag are dependent on 


for a double parabolic arc air 
the thickness ratio, hypersonic similarity parameter and normalized 
angle of attack. A chordwise and a normal force coefficient are 
generated. Lift and drag are found from these values. Also included in 
drag is skin friction drag calculated for the wing and body, as 
previously stated in the subsonic and supersonic sections. The values 
agree well with experiment for wings of finite aspect ratio at small 
angles of attack. This is because the region of three-dimensional flow 
occurs only over a small portion of the wing planform at hypersonic 
speeds ("Hypersonic Aerodynamics" p.63). 

PROGRAM USE 

For the final configuration 
shown in figure 3.1, the aircraft was 
modeled as a delta wing with a small 
body in the center for fuel. This 
area is shaded in figure 3.2. The 
areas and lengths input to the pro- 
grams are given in figure A. 9. Other 
input data included an average thick- 
ness to chord ratio of 3/39=9.1, a 
leading edge radius of .92% (which equal a 1/8 inch radius, which 
according to temperatures estimation will be sufficient), a sweep angle 
of 76 degrees and body width, and body cross-sectional area. 

Using this data in the programs, the pertinent aerodynamic was 
obtained. The non-dimensional representation of this data follows this 






sec t i on . This data was then used by 
the propulsion group along with the 
engine data to obtain a viable flight 
profile. The glide portion of the 
flight profile was determined by the 
aerodynamics group since no engines 
were used in this portion. A brief 
description of this calculation will be 
given • 
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Aerodynamic Results 


Cl vs Angle of Attack 


/ 

/ 

_ 

/ 

/ 

/ 

M = 0.3 j 

/ 

M = 2.0 | 

/ 

t 

s 

i 

1 M = 5.0 | 

s 

s 




// 


/ 

/ ' 


' 


6 0.5 i L5 2 25 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 3 3.5 9 9.510 

(X (deg) 



C L vs Mach Number 



ORiGiNAL GAGE IS 
OF POOR QUALITY 











Chapter 4 

LONGITUDINAL STABILITY 
(John Aslan) 

The longitudinal stability of the aircraft must be known throughout 
the entire flight to determine what type of, if any, stability augmenta- 
tion system must be implemented. The reference used is "Fundamentals of 
Aircraft Design" by Leland M. Nicolai, 1975. The basic method of this 
reference is used in the analysis, with some new assumptions to tailor 
it to that of a wing-body design. The contribution of the wing-body and 
propulsion units to the moment about the center of gravity will be used 
to determine the static longitudinal stability. 

The two sets of winglets at the rear of the aircraft were found to 
contribute about 2/3 of the lift of the aircraft, and the rest of the 
aircraft as 1/3 of the lift. This was determined by finding the lift- 
curve slopes of the individual parts by the same method that was used in 
determining the lift-curve slope of the entire aircraft (see Aerodynam- 
ics) . The aerodynamic center of the entire aircraft was then assumed to 
be that of the winglets. The aerodynamic center of the entire aircraft 
and the rest of the aircraft were found to vary between 40% and 50% of 
their respective mean aerodynamic chords, depending on the flight mach 
number . 

Now that the aerodynamic centers of the aircraft have been found, 
the center of gravity and its travel must be found. This was accom- 
plished using a weight estimation program, which is described in this 
report. Next, a basic moment equation is written, taking into consider- 
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at ion the thrust vectors and lift and drag at arbitrary angles of 
attack. The longitudinal distances are referenced to a fixed point on 
the aircraft, namely the leading edge of the winglets. To simplify the 
equations, the following assumptions, were made: 

1. small angle of attack 

2. vertical contributions negligible, because they are 
relatively small distances compared to the horizontal 
di stances 

3. lift coefficient >> drag coefficient 

4. the moments about the aerodynamic centers of the 
aircraft components are not contributing, because no 
suitable method for determining these was found. 

From the final moment equation, the moment coefficient at zero 
angle of attack and the moment-curve slope are found by breaking down 
the single moment equation into a constant part and one that depends on 
angle of attack. Once these equations are known, the center of gravity 
and aerodynamic centers are determined for a specific time in the 
flight. The distances and forces are substituted into these equations, 
and the curve found. 

During the turbo fanramje t only phase of the flight, landing (plot 
at end of chapter) and acceleration from Mach 2 to Mach 5. the aircraft 
is statically stable. However, with the scramjets or rockets running, 
the nose-down mo men t is quite considerable, which makes the mo me n t at 
zero angle of attack negative, with the slope still negative. In this 
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configuration, trimming is not possible. A stability augmentation 
system will be needed to accomplish this. Also, a plot at the end of 
the chapter shows the variation in static margin as the fuel is used in 
tank order A, C, then B. As can be seen, the static margin remains 
almost constant and is positive, which is a criterion for static 
longitudinal stability. Further analysis and wind tunnel testing are 
needed to get an accurate assessment of the situation. 
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Chapter 5 


SONIC BOOM 
(Terry Bisard) 

Int roduc t i on 

With the mission profile requiring high supersonic flight, the 
amount of sonic boom produced by the aircraft becomes very important. 
A sonic boom is produced when the plane travels at supersonic speed and 
forms shock waves. A large rise and fall in pressure results from the 
quick passage of the shock waves. These sonic booms passing over land 
can be strong enough to damage people’s hearing and even shatter 
building windows. Due to this, the F.A.A. has stipulated that no 
aircraft may produce a sonic boom over 1 lb/ft~2 while traveling over 
United States airspace. 

Method of Calculation 

In designing our supersonic plane, flight path consideration must 
be given to the sonic boom it produces. The Morris prediction method 
was used to estimate the sonic boom signature produced by our aircraft. 
This method assumes a flying diamond shape and the change in pressure is 
given by: 

dP 1 = . 363*Kr* ( Pa/Pg) ~ . 5*K1 * ( M~2- 1 ) ~ . 375/M*W~ . 5*Lw~- . 25*Pa~ - . 5* 

H~- . 75*Pg 

Where: Kr , Kv - Shape factor constants 

Pg,Pa - Pressure at altitude and sea level 
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M - Mach number 

W - Weight of aircraft 

Lw - Length of aircraft 

H - Flight altitude 

Results 

By writing a computer program, the sonic boom signature can be 
found. Using this method, the sonic boom can be found for the flight 
profile at various alti- 
tudes and Mach numbers. 

As is seen in fig. 5.1, 
our sonic boom signa- 
ture was very small and 
was above legal limits 
only for Mach numbers 
between 1.1 and 2.5. 

Even when our sonic boom 
was above the legal lim- 
it it did not fan out 
very far. The small 
sonic boom signature 
comes about from the 

, , , . Figure 5.1: Surface Boom Pressure 

aircraft’s relatively 1 

low weight ( ~5Q , 000 lb.) and a short total length of only 60 feet. 




Chapter 6 


PROPULSION 
(Warren Peters) 


I nt roduc t i on 

The ultimate purpose of the propulsion system of this aircraft is 
to accelerate to the desired mach number and to test the scramjet 
engines. One of our personal goals was to keep the plane as small as 
possible, which means a compact power plant configuration and the 
minimum amount of fuel necessary to complete the mission. 

Of the many types of propulsion systems available, we decided to 
use a combination of turbof an/ramjets , scramjets , and rocket engines. 
Liquid rocket engines were eliminated from consideration because of the 
large weight of the hardware necessary to use one, so the propellant of 
choice was a solid propellant. Calculations show that solid rockets are 
too heavy and bulky to use on our small aircraft if it is the only 
propulsive system available, as shown last quarter. However, if the 
solid rockets are used to augment a turbof an/ramje t , the results are 
quite acceptable. The turbofan/ramjet engine is capable of producing 
most of the thrust necessary to accelerate from launch to M = 6, provide 
a powered landing and is insurance in the case of a scramjet failure. 
The scramjets are necessary for the plane because testing them is the 
purpose of this analysis. 

Last quarter we ended with a propulsive configuration of 2 
scramjets and 1 turbofan/ramjet engine that was sized down to 90%. The 
total amount of fuel used in this configuration was 12,000 lbs, which 
was far too large for the size of the plane we wanted this quarter. The 
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total amount of fuel contained in this quarter’s design was 7,100 lbs, 
and the power plant system was designed around this. An extra scramjet 
engine was added this quarter to bring a total of 3 scramjets. More 
scramjet engines were not used because the small size of our plane would 
not permit it. Two small external solid rockets were also added to the 
design that assist in the thrust pinch and low speed operations, and are 
ejected after they have burned out. The third and final change was that 
the t urbo f an/ r am j e t was upscaled by 2%, leaving the total dimensions at 
92% of full size. 

Rocket 

The largest innovation in the powerplant package was the addition 
of 2 solid rocket engines that are used in the early stages of the 
flight. As mentioned earlier, solid rockets were dismissed for use as 
the only propulsion system from lanuch to M = 6 because of the severe 
weight penalty incurred. However, a small set of solid rockets used in 
conjunction with the turbofan/ramjet engine increases performance and 
decreases the hydrogen fuel weight. The following is a list of the 
major advantages of using the solid rockets. 

Performance - The high thrust levels produced by a solid rocket 
engine increases the acceleration loads to values greater than lg. This 
large acceleration decreases the time that the aircraft burns fuel from 
launch to M = 6. Since the turbofan/ramjet engine burns less hydrogen 
fuel for a smaller amount of time, the total amount of hydrogen fuel 
used decreases. The decrease in hydrogen is necessary because of it’s 
low density and large volume requirements. Another problem is that the 


calculated transonic drag values come very close to the levels of thrust 
available from the t urbo f an/ ram je t at the thrust pinch. The innovative 
configuration of our aircraft makes it difficult to accurately describe 
the drag levels, but if rocket engines are used, the thrust levels are 
extensive. This increase in thrust is much more that the increase in 
drag, and the worry of the thrust pinch has been put to rest. 

Size/Weight - The rockets used on our aircraft are quite compact 
compared to the airbreathing systems and fuel. The density of the solid 
fuel makes the rockets small enough that the drag is not an inhibiting 
factor. Since the rockets are attached to the exterior of the aircraft, 
drag is a consideration. Although the drag is by no means negligible, 
the increase in thrust is more than enough to compensate. Another 
attractive feature of placing the rockets on the exterior is that the 
plane can remain small since it does not need to use the volume on the 
interior of the aircraft. The weight of the casing and components 
necessary to contain a solid rocket are small when compared to either 
the we i gh t of solid propel 1 an t or the hardware that would be utilized to 
ran a liquid rocket engine. 

Expendabi 1 i ty - The result of placing the rockets on the exterior 
of the aircraft and using solid propellants is that the rockets may be 
ejected after they have burned out. The rockets will produce drag as 
thev are being used and will add to t he total mass of the aircraft, but 
may be ejected after their purpose has been fulfilled. The drag 
producing shells and extra weight of the hardware will be absent from 
the larger mach number flight. 
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Design 


A particular solid rocket that would fit our specifications is not 
already in production, so a rocket was designed to fit our needs. The 
solid fuel that was used is the same fuel that is used in the first 
stage of the Minuteman Missile. The fuel specs were obtained in Rocket 
Propulsion Elements by George Sutton, see references. The following 
table is the list of the fuel specifications and the final design 
di mens ions . 


Fuel composition 16% Al, 70% NH4CI04 


Dens i ty 
Burning Rate 
Specific Impulse 
Outer Diameter 
Length 
L/D 


14% polymer binding 
.0636 lb/inches cubed 
.349 in/s 
214 s 
30 in 
12.33 ft 
4 . 9 


Grain Configuration - Slotted Tube 


Chamber Pressure 850 psia 
Burn time 29.8 s 

Thrust (constant) 25000 lbs 

Fuel Weight 3480 lbs 

Web thickness 10.4 in 

Elastomeric liner .1 in 
case weight 113 lbs 


Rocket Design Calculations 

The calculations that were used can be seen in figure 6.2. These 
calculations were based on the process found in Rocket Propulsion 
Elements. The criteria set out from the aircraft is that one rocket 
must produce a constant 25000 lbs of thrust for 29.8 seconds. Also, the 
length must be less than 13 ft in order to fit in the correct place on 
t he aircraft. 
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Figure 6.2: Rocket Equations 


Rocket Spreadsheet Design 

The equations as set forth in figure 6.2 have been entered into a 
spreadsheet program. See figure 6.3. I varied the parameters of total 
thrust and outer diameter, and the outputs were weight and length. 
Length/Diameter was also calculated so that the correct grain configura- 
tion could be determined. The final dimensions of the rocket have 
already been displayed in the table contained in the Rocket Design 
Section. 
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Rocket Design Program 


Force 

Ou t s i de 
Diameter 

Burn 

Time 

Weight 

Vo lume 

Ins i de 
Diameter 

Length 

L/D 

( lbs ) 

(in) 

(s) 

( lbs ) 

( in) 

( in) 

( in) 


30000 

50 

49.68944 

3387 .916 

54334.51 

15.31677 

53 . 33614 

1 .066723 

30000 

48 

49.68944 

3387.916 

54334.51 

13.31677 

62.83041 

1 . 308967 

30000 

46 

49.68944 

3387.916 

54334.51 

1 1 .31677 

75 . 76706 

1 .6471 1 

30000 

44 

49 . 68944 

3387.916 

54334.51 

9.31677 

94 .37102 

2 . 144796 

30000 

42 

49 .68944 

3387.916 

54334.51 

7.31677 

123 . 3009 

2 .935737 

30000 

40 

49.68944 

3387.916 

54334.51 

5 .31677 

174 . 2269 

4 . 355673 

40000 

44 

37 .26708 

3387.916 

54334.51 

17.98758 

54 .93365 

1 . 248492 

4Q0QQ 

42 

37.26708 

3387.916 

54334 .51 

15 . 98758 

63.62317 

1 .514837 

40000 

40 

37.26708 

3387 .916 

54334.51 

13.98758 

74 . 9235 

1 .873087 

4 GQQQ 

38 

37.26708 

3387 .916 

54334 .51 

11 .98758 

90 . 15515 

2 . 372504 

40000 

36 

37 . 26708 

3387.916 

54334.51 

9.987578 

1 1 1 .6985 

3 . 102736 

40000 

34 

37 . 26708 

3387.916 

54334.51 

7.987578 

144 .3212 

4 . 24474 

40000 

32 

37 . 26708 

3387.916 

54334.51 

5.987578 

199 . 1655 

6 . 223921 

5QQQQ 

44 

29.81366 

3387.916 

54334.51 

23. 19006 

46.03018 

1 .046141 

50000 

42 

29.81366 

3387 .916 

54334.51 

21 .19006 

51 .97873 

1 . 237589 

50000 

40 

29.81366 

3387.916 

54334.51 

19 . 19006 

59 . 2837 

1 . 482092 

50000 

38 

29.81366 

3387.916 

54334.51 

17 . 19006 

68.43182 

1 .800837 

50000 

36 

29.81 366 

3387.916 

54334.51 

15 . 19006 

80 . 16824 

2 . 226896 

50000 

34 

29.81366 

3387.916 

54334.51 

13. 19006 

95.69298 

2 . 8 14499 

50000 

32 

29.81366 

3387.916 

54334.51 

1 1 . 19006 

1 17.068 

3 .658374 

50000 

30 

29 .81366 

3387.916 

54334.51 

9. 190062 

148.156 

4 .938532 

50000 

28 

29.81366 

3387 .916 

54334.51 

7 . 190062 

197 . 1266 

7.040237 

60000 

40 

24 .84472 

3387.916 

54334.51 

22 .65839 

53 . 24605 

1 . 331 151 

60000 

38 

24 . 84472 

3387 .916 

54334.51 

20 . 65839 

60.51 152 

1 .592408 

60000 

36 

24.84472 

3387.916 

54334 . 5 1 

18 . 65839 

69 . 5098 

1 .930828 

60000 

34 

24.84472 

3387.916 

54334.51 

16.65839 

80 . 88795 

2 . 379057 

60000 

32 

24.84472 

3387 .916 

54334.51 

14.65839 

95 .65039 

2 . 989075 

60000 

30 

24.84472 

3387.916 

54334.51 

12.65839 

115.4423 

3 .848076 

60000 

28 

24 .84472 

3387 .916 

54334.51 

10 .65839 

143.1521 

5.112576 

60000 

26 

24.84472 

3387 .916 

54334.51 

8 .658385 

184 . 3505 

7 . 090403 

60000 

24 

24 .84472 

3387 .9 16 

54334 . 5 1 

6 .658385 

251 .3216 

10.47173 


Figure 6.3: Rocket Calculations 


Turbofan/Ramjet (M=.8 to M=6, and landing) 

For the acceleration from launch to M=6, 1 t u rbo f an/ ram j e t engine 
sized to 92% of full scale was used. The reason for scaling the engine 
down was to integrate it correctly into our airframe configuration. 
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There were two main purposes in using a turbofan ramjet in our 
configuration. The first reason is that we wanted a powered landing. 
In the case that the calculated flight trajectory would not match the 
test conditions, power may be needed to maneuver the aircraft at or near 
the time of landing. Another prominent reason is that if the scramjets 
should fail at anytime during the test, the turbof an/ ramje t would be 
able to provide enough power to get the aircraft to Edwards. There is 
enough fuel for the turbof an/ramjet to take the plane the entire 
distance from launch and still have excess fuel left over. As a final 
note, since no other propulsion system could be found that provided data 
for analysis, we has little choice but to use the turbofan/ramjet. 

Data 

Detailed analysis of the performance of the turbof an/ramj e t engine 
was given to us by General Electric. The scaling relationships for the 
engine were also given to us by G.E. The data consisted of thrusts, 
fuel flow rates ( FFR ) and airflow rates at different altitudes, mach 
numbers and SFCs . The data was not readily analyzable until it was 
transformed into a more usable form. The data was somewhat sparse, so 
I developed a scheme to relate these points into a continuous curve that 
could be used to find thrusts and FFRs at any altitude and mach number. 

The first step was to plot the thrust versus altitude on a graph 
using the few points given to us. A curve was then fit through these 
points to provide a continuous thrus t /a 1 t i tude relationship. This was 
plotted for 7 different mach numbers that would be encountered through 
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the operating region of the 
turbofan/ramjet operation 
(M=.8to6). See figure 6.4. 

The next step was to 
repeat this process on anoth- 
er graph , but to use FFR in- 
stead of thrust in order to 
obtain a FFR vs thrust curve. 
Again, this was given at dif- 
ferent mach numbers as shown 
in figure 6.5. Both the 
thrust and FFR calculations 
were given for the maximum 
attainable SFC at each point. 
This was chosen because anal- 
ysis showed that the larger 
acceleration from this higher 
thrust overcame the increased 
FFR and gave the minimum 
amount of fuel burned. 

The final step was to 
use the values from the aero- 
dynamic program that gave 
drag at different altitudes 
and mach numbers. This set 
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Figure 6*5: Fuel Flow Rate 


of data, drag versus Altitude at a given mach number, was overlaid onto 




the thrust vs. alt graph because both of these forces are used in con- 
junction in the calculation of the flight profile. See figure 6.4. 

Scramjet (M=6 to M=9.25) 

The purpose of this vehicle is to test a scramjet engine at a high 
mach number after accelerating from M = 6 at 100k ft. The total number of 
scramjets chosen was 3. The size of the scramjet could not be changed 
since no scaling relationships could be found. Instead of scaling the 
engine, multiple engines may be used as a single propulsive unit. We 
used 3 scramjet engines because this was the size limitation of our 
aircraft. Any more engines would have improved the fuel rates, but 
would have decreased the lifting surface and destroyed the aerodynamics 
of the upper lifting surface. 

Using the scramjet data that was provided by GE for their engines, 

I again set up a set of curves that would be the relationships for a 
more continuous distribution of thrusts vs speeds. I mention speed only 
because the scramjets will operate only at lQQk ft, thus, the altitude 
is not a variable. I used the GE curves of thrust vs fuel equivalence 
ratios ( ER7 ) at different speeds and replotted it as thrust vs speed at 
a different ER7 values. See figure 6. This was preferable because it 
provided a thrust output for varying speeds, which are the parameters 
thcit the fuel program uses as inputs. Please note that drag was a given 
distribution over the M=6 to M=1Q speed range due to constant altitude. 
The drag was not plotted on the thrust curve because of the simple 
nature of the trajectory. 
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The fuel flow rate vs speed was also plotted for the scramjet 
engines using the method mentioned for the thrust levels. See figure 

6.7. For each ER7 , there was a 
given thrust and FFR, so as can be 
seen, the ER7 is the unifying fac- 
tor between the graphs. 

Turbofan/Ramjet-Scramjet Fuel Pro- 
gram 


The purpose of this program was 
to take the thrust, drag and fuel flow 
rates and calculate the amount of fuel 
used, distance, accelerations etc.. 
This program was used for the turbo- 
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Figure 6.7: Scramjet Fuel Flow 
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Figure 6.6: Raw Scramjet Thrust Data 


fan/ramjet and scramjet acceleration 

plus the scramjet cruise time. The glide portion and subsonic cruise 
before landing are handled in the glide portion of this report. The 
program is based on Newtons Second Law and the relationship between FFR 
and time. The program was integrated into a spreadsheet so that the 
inputs could be easily changed and the results readily graphed. An 
example of a portion of the numerical spreadsheet listing is located in 
Figure 6.8. A few of the particulars are included and list as follows: 
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Scaled 90% 


l Engine 44000 


Ml Total Total FFR 
Thrust lbs/s 


0.9 

66928 

3.30096 

1 

74122.4 

4.57056 

1 .5 

82586 . 4 

6.47496 

2 

89780.8 

8 . 2524 

2.5 

46552 

9.9452 

3 

49937 . 6 

1 1.638 

3.5 

52900 

13.9656 

4 

52900 

14.812 

4.5 

48244.8 

13.5424 

5 

51111.2 

16.8948 

5.5 

37568.8 

13 . 5092 


Scrarajet Acceleration 
Ml 


6 

7 

8 
9 


Scramjet Steady State 
M 


9 .25 


Drag 

Acc 

Time 

Fue 1 

Weight 

lbs 

ft/s 

s 

lbs 

1 bs 

19000 

35.89 

2.7027 

8.9214 

4359 1 

21000 

41.72 

1 1.625 

53.133 

41338 

23000 

49.2 

9.8583 

63.832 

39074 

27000 

53 . 2 

9.1168 

75 . 236 

36799 

21500 

21.92 

22.125 

220.03 

36579 

23900 

22.92 

21.16 

246 . 26 

36333 

25500 

24 . 28 

19.972 

278 . 93 

36054 

25000 

24.92 

19.464 

288 . 3 

35765 

23000 

22.73 

21 .339 

288 . 98 

35476 

18500 

29.6 

16.385 

276.83 

35200 

11270 

24.06 

20.16 

173.91 

272 . 34 
2072 . 8 

34927 


Drag 

ACC 

Time 

Fue 1 

We i gh t 

lbs 

ft/s 

s 

lbs 

lbs 

8500 

5.26 

190 . 1 1 

1209 . 1 

33990 

9500 

6.915 

144.6 

1214.7 

32776 

10400 

8.94 

111.86 

1211.4 

3 1564 

11001 

11.42 

21 .883 

296 . 07 

31268 



468 .45 

3931 . 3 



Drag ACC Time Fuel Weight 

lbs ft/s s lbs lbs 

1 1 300 0 120 7 20 30548 


Fue 1 
1 bs 

6724 . 05 


Figure 6.8: Turbofan/Ramjet Spreadsheet 
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- The program is based on a changing mass as the fuel is continu- 
ously being burned. 

- Using Newtons Second Law, the thrust, drag and mass are used to 
calculate accelerations, distances and times derived from 

the kinematic equations of motion. 

- Uses changing FFRs and time in order to calculate the fuel usage 


The need for this pro- 
gram and the versatility 
built into it will now be 
discussed. The beginning of 
our iteration loop for the 
turbofan/ramjet starts with 
the designation of a flight 
profile (M vs Alt). The 
thrusts and drags can now be 
obtained from the raw data 
curves in figure 6.4. Using 
the same flight profile, the fuel flow rates can be obtained from the 
graph in figure 6.5. The thrusts, drags and fuel flow rates can now be 
easily entered into the program. The outputs are distance, time, fuel 
and accelerations. Since the fuel weight is our major concern, we 
continued to pick various flight profiles until we optimized the fuel 
burned in this region. The difference between the best and worst flight 
profiles was approximately 4000 lbs, which is half of the entire fuel 
weight! The optimization of the profile was critical in the success of 
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Figure 6.9: Altitude vs. Mach No. 
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our aircraft. The final result of this profile iteration is shown in 
figure 6.9. 

The iteration loop for the scramjets was somewhat different. A 
flight profile was not the beginning of the iteration because we were at 
a constant altitude and the profile was trivial. The start of this 
iteration loop was the choice of a single ER7 or a combination of 
different ER7s over the speed range used to read thrusts and FFR from 
the graphs in Figure 6.6 and Figure 6.7. These values were fed into the 
scramjet portion of the program Again, the fuel weight burned was the 
prime consideration. The result of this optimization was that a 
constant ER7 of 1.2 was to be used for the entire speed range. This 
saved us another 20QQ lbs of fuel over the case of a poorly optimized 
cho ice of ER7 . 

The scramjet cruise portion was included on the spreadsheet but 
could have been done by hand. The thrust was set equal to drag and the 
corresponding FFR was used for the 2 minute cruise. 

Profile (Acceleration and Cruise) 

All of the graphs that will be presented in this section are direct 
results from the fuel spreadsheet program. All of these cases were 
considered for both the configuration that included the rocket and that 
which was ’’rocket less" . 

Fuel Percentages 

See figures 6.10 & 6.11. The comparison of the rocket and 
rocket less versions both show that over 50% of the fuel is used in the 
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scramjet acceleration. The turbo- 
fan/ramjet cruise is the other 
major portion of fuel. Note that 
adding the rockets decreases the 
amount of fuel needed in the tur- 
bof an/ramjet and lets the scram- 
jets do more work and attain a 
higher mach number. 


FUEL PERCENTAGES 

Internal Fuel = 7100 lb 
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Figure 6.10: Rocketed Fuel %. 
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FUEL PERCENTAGES Force vs M 

Internal Fuel = 7100 lb See figures 6.12 & 6.13. This 

graph plots the thrust and drag vs M 

• for the entire acceleration stage. 

The primary design is to keep the 
^ 2750 ^ thrust curve higher than the drag 
curve at all times. On the graph of 
the rocketless configuration in fig- 
ure 6.13, you can see that at the 
B transonic region the aircraft ex peri- 

ences a thrust pinch. This is due to the high drag because of transonic 
effects on aerodynamics. The rocket version in figure 6.12 shows our 
solution to this problem. The high thrust is considerable compared to 
the smaller increase in drag. The result is not only less fuel used, 
but there is no fear of not breaking the speed of sound. 


(3270)-^||P^ 

Figure 6.11: Rocketless Fuel % 




Acc vs M 


See figures 6.14 & 6.15. These 
graphs depict the data found in the 
spreadsheet program. The rocket 
accelerations are considerably higher 
in the earlier region as would be 
expected. The negative accelerations 
are due to the glide portion of our 
profile (deceleration). The slight 
discontinuity near M=5 to 6 is due to 

the fact that the scramjets were 
turned on a little early to aid in 
the acceleration at that point. 

The rest of the figures (6.16 
through 6.21) are plots of different 
variables through the profile and 
have been provided to give any addi- 
tional information that you may have 
wanted to see . 



Figure 6.14: Rocketed Acc vs. M 


FORCE vs. MACH NUMBER 

Rocket Configuration 



Figure 6.12! Rocketed F vs. M 
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DISTANCE vs. TIME 
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Figure 6.21: Dist. vs. t 
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Chapter 7 


GLIDE CALCULATION 
(Gary Huf ford} 

The last portion of the mission profile is the glide and landing 
portion. After the two minute Mach IS cruise has been completed, all 
engine systems are shut down. The airplane will then gradually decrease 
Mach number and altitude until Mach Q.9 at 4Q,SGQ feet is reached. This 
point will also be thirty miles off the California coast in preparation 
for landing at Edward’s Air Force Base. This point was chosen such that 

the sonic boom created as the airplane reaches Mach Q.9 from Mach IS 

will not affect the California coastline. 

After this 75 mile cruise the airplane will begin its landing 
descent. This descent will be gradual from Mach S.9 at 4Q,QQS feet to 
Mach 0.3 at sea level. Thirty miles have been allowed for this descent. 
This last 3Q miles will allow landing at Edwards Air Force base. 

' The fuel needed for the Mach Q.9 cruise is 227 pounds. This was 

computed by setting thrust equal drag and finding the corresponding fuel 
flow rate. By computing the velocity, the time needed for a 75 mile 
cruise was calculated. The fuel was computed by multiplying time and 
f ue 1 f l ow rate. 

The fuel needed for the thirty mile descent to landing is 132 

pounds. To compute this value the deceleration to travel thirty miles 
must first be computed. From the acceleration, time and thrust were 
computed. With the thrust level known the fuel flow rate is determined. 
Fuel needed is then computed by multiplying fuel flow rate and time. 
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The glide descent and fuel computation are given in tabular form in 
Chapter 6. The glide and landing profile were chosen arbitrarily to 
obtain approximately a 5QQ mile distance. This portion of the mission 
profile can be seen graphically in Chapter 8. The fuel computations are 
given in Chapter 6. 
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Chapter 8 


FLIGHT PLAN 
(Gary Huf ford) 

The flight plan is shown in Figure 8.1. The carrier aircraft will 
take-off from McChord Air Force Base in Washington and fly over the 
Pacific Ocean where the vehicle will be launched. The research vehicle 
will fly in a straight line through the Mach 10 cruise and descent until 
the aircraft reaches Mach .9 at 40,000 feet approximately 30 miles off 
the coast. At Mach 0.9 and 40,000 feet approximately thirty miles off 
the coast. At Mach 0.9 and 40,000 feet the subsonic cruise will begin 
to land 'at Edwards. 

Four reasons exist for choosing this flight plan. First, we must 
fly over water for safety reasons. If something should go wrong, no 
lives will be endangered. Second, the length of the flight plan is the 
projected range of the vehicle for a successful test. At any point 
during the test, if the test is aborted sufficient fuel will be left to 
land at Edwards Air Force Base correctly utilizing the turbo fanramje t . 
Third, we must fly close to land to be able to maintain ground based 
communications. At Mach IQ the vehicle will cover large distances and 
several ground based stations may be necessary to maintain constant 
communication. And fourth, the sonic boom problem becomes non-existent 
for this flight plan. By flying through Mach 1 at 40,990 feet, thirty 
miles off the coast, the sonic boom produced is well within tolerable 
limits. 
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Figure 8.1: Flight Plan 



Chapter 9 


SCRAM JET INLET DESIGN 
(Terry Bisard) 

Introduction 

The purpose of the inlet is to slow down the free-stream flow and 
to create a favorable pressure gradient and weight flow rate to enter 
into the combustion chamber. Designing a supersonic inlet is a very 
complicated task. In actuality an inlet is of a three-dimensional 
design, but to make caculations simpler, our inlet assumes a two 
dimensional wedge shaped design. Using oblique shock theory and the 
method of charac te i s t i cs , the flowfield inside the inlet can be 
ana 1 yzed . 

For our supersonic research vehicle, we decided to put our 
SCRAMjets on the top surface of the aircraft and the turbo fanramje t on 
the bottom of the aircraft. This allows two separate inlet systems, to 
be designed to make the two different engine systems to operate more 
efficiently. Also, because the supersoic inlet only has to vary from 
Mach six to ten, the total inlet geometry is less complex. Inlet 
Geome t ry 

The compression angle is the ang'le at which the inlet converges* 
Values for the compression angle vary between six and eight degrees. A 
few drawbacks to high compression angles are high total pressure losses, 
and early shock induced boundary layer separation. The only problem 
with a low angle is increased inlet length. As the compression angle 
increases the flow gradients become steeper. This occurs because the 



of these conditions must be met while keeping the inlet length less than 
18 feet. 

Flowfield Analyzation 

In analyzing the inlet, one must consider the series of reflected 
shocks that develop. Because of these shocks, the flow is slowed down 
and there is a tremendous loss of total pressure. Instead of lengthy, 
complex hand calculations, a computer program was developed to calculate 
the flowfield parameters at certain stations along the inlet. Oblique 
shock analysis can be used because our SCRAMjet inlet assumes a two- 
dimensional wedge shape design. 

To maintain organization, the main program consists of many program 
subroutines. The OBLIQUE subroutine calculates flow conditions behind 
an oblique shock wave and the GEOM subroutine pinpoints the position in 
the inlet that these conditions occur. The program should be run for 
several different cases to find the best inlet design configuration. 

Boundary Layer Interaction 

In analyzing any body in a hypersonic flow, the viscous boundary 
layer must be taken into consideration. This is especially important in 
inlet design, because turbulent flow entering a conbustion chamber could 
disrupt an efficient thrust process. For a SCRAMjet inlet, the boundary 
layer is very large and is on the order of 20% or more of the entile 
flow entering the combustor. 

To insure that only a small boundary layer interacts with the flow 
■foin^ inot the conbustor, a boundary layer reduction system needs to be 
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Chapter 10 


TURBOFAN/RAMJET INLET 
(Jin Sogers) 

Unlike the Scramjet Inlet which operates only in a narrow mach 
number region, the turbo fan/ ramje t inlet must operate from subsonic 
speeds from about mach 0.4 during landing up to a blistering mach 6.0. 
The accomplishment of such a design is not a trivial task. The 
requirements of the system change drastically throughout the operational 

region. 

Through each region, the inlet must provide the appropriate mass 
flow required by the engine at the correct conditions. Along with this, 
the air flow must be of good quality. That is to say, it should have 
approximately the same conditions across the compressor. If inequali- 
ties exist, the results could be disastrous. 

DESIGN 

The first concept desired for this aircraft was a conical inlet. 
This was mainly for aesthetics, but never- the- 1 ess it was the first. It 
soon became apparent that a conical inlet would not suit our purposes. 
In order to accomplish the mission throughout the profile, the inlet 
would have been far too complex to be feasible. The next logical step 
was to design a two-dimensional inlet. 

By studying the required ramp angles to provide a suitable pressure 
recovery at mach 6, a few tendencies were noted. One of these was that 
if the ramps generally increased as the mach number decreased through 
t}-e channel, the pressure loss wouldn’t be quite as bail. At this point 
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as much as 2 inches at the inlet. , 

Due to the unknown nature of the | — ~~ ^ 

boundary layer in the hypersonic 

28 ’ 

region, a 4 inch initial bleed r^. \ ^ 

diverter is implemented just prior 1 • ^ 

to the inlet. If the boundary / 

' .74 

layer becomes too much of a prob- Fiffure 10 . 2: External Ramps 

lem in this area, surface cooling 

will be app lied in order to shrink it. This doesn’t eliminate the 

boundary effects. The boundary layer must be stripped throughout the 

inlet. The external ramp section will be made of perforated metal 

sheets which can provide just the right amount of mass flow to strip the 

boundary layer as it attempts to grow through this section. To account 

for all boundary layer mass losses and possible auxiliary system losses, 

2 

the actual capture area is increased by 2®%, to 6 ft , 



smooth increase in area at all mach numbers. The 2-D section of the 
inlet is shown in figure 10.4. 

At the end of the 2-D section of the inlet, another system of 
boundary layer diverters is located. At this point a rectangular to 
circular divergent section will gate the flow into the compressor. This 
section is one solid piece, except for two 4 ft door sections on either 
side which provide the necessary mass flow rate at the subsonic speeds. 
The engine will be in a starvation period through some of the faster 
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Figure 10.5! Turbofan/Ramjet Inlet 3-view 



64 


: !LMED 



subsonic speeds. This is acceptable due to the rocket assist on ascent 
and the qu ick descent, where this is though t not to be critical. 

The test of efficiency for an inlet seems to be the level of 
pressure recovery available. Figure 19.6 shows the pressure recovery 
throughout the flight profile. 






Chapter 1 1 


COOLING 

(Jim Rogers) 

Int roduc t ion 

As a body moves at low supersonic speeds, the heating of the body 
is not a very important consideration. As the speed increases, the 
heating of the body increases. At hypersonic speeds, as those encoun- 
tered by the vehicle under study, the heating becomes a primary concern. 
For example, the recovery temperature at Mach 10 is greater than 4Q0Q°R. 
There are very few materials which are solid at these temperatures, and 
none (as of yet) that can be used reliably for an aircraft’s structure. 

Even if these materials were available, the vehicles internal 
temperature would be far too great for the electronics, fuel system, and 
hydraulic system. If the vehicle were to be manned, a further con- 
straint on the internal heat would be present. The pilot would also not 
be able to exit the aircraft until after a lengthy cool-down period. 
This design, being unmanned, will not be concerned with post-landing 
problems at this time. 

Cooling Methods 

As a result of the previous discussion, this aircraft will need to 
have some sort of cooling system. There are two basic types of cooling, 
active and passive. Passive cooling attempts to design the structure in 
such a way so as to provide a highly conductive path to shunt the 
heating away from the "hot spots" to the cooler areas where the heat may 
be radiated away. This type of system cannot provide the cooling rate 
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required near the forward surfaces, therefore, although some passive 
cooling may be found to be useful, the subject will not be discussed 

further here. 

Since passive cooling is not adequate, an active cooling system 
needs to be considered. By using some secondary fluid, the heating can 
be diverted into this fluid and carried away. There are three basic 
types: transpiration, direct convection, and indirect convection. 

Figures 11.1 thru 11.3 display schematics of each of these systems. 



Figure 11.1: Transpiration cooling 



Transpiration cooling, illustrated in figure 11.1, uses a porous 
surface material with which a fluid is forced to seep out to the 
boundary layer. This provides a layer of fluid in between the hot 
boundary layer and the cooler aircraft surface. The heat is still 
transferred from the boundary layer, but the surface into which it is 
going is being continuously replaced. This mass transfer into the 
boundary layer causes the boundary layer to increase in size while 
decreasing the surface friction. This system appears very good 
aerodynamical ly, but does require the aircraft to carry the fluid which 
increases overall weight. This fluid could be water (actually some 
water mixture such as water-glycol), hydrogen or just air. 

Air is appealing because it does not need to be carried in the 
aircraft. The air would be "scooped up" as the aircraft is flying. 
This method may naturally provide the necessary pressure to drive the 
s vs t em . The primary problem with this method is the capture area 
required to provide the necessary mass flow rate. If the mass flow can 
be provided, the air can be passed through a heat exchanger to cool and 
then passed through the porous surface providing the cooling. The heat 
exchanger cools the air by heating the hydrogen fuel on its way through 
the fuel system. 

A variation on this theme is to use water or hydrogen to provide 
the fluid cooling layer. The advantages of these methods is that the 
fluid does not need to be cooled before used and the lack of the large 
capture area reduces drag. The disadvantage is the additional weight of 
these fuels. Although hydrogen may be lighter, the increase in volume 
and the safety issues preclude its use. Another problem is the 
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streaking which may occur, possibly leaving the surface dangerously 
exposed to the extreme conditions it wasn't designed to withstand. 



Figure 11.2 illustrates the direct convection method. This method 
directs the hydrogen fuel around the aircraft cooling the surfaces by 
transferring the heat to the fuel. This method is an effective cooling 
technique but does have its problems. The dangers of an aircraft 
literally lined with a dangerously explosive substance needs no 
elaboration. Another problem is hydrogen's tendency to react with 
titanium. Because the structure, at this point of the development, will 
be made of a titanium-aluminum alloy, this method looks unattractive. 

A variation on this method is called spray cooling. A pressurized 
coolant is sprayed through a nozzle causing the liquid droplets to 
vaporize as the hit the hot surface. The resultant vapor is routed to 
another portion to be cooled and then pumped back for repeated use. 
This method can use a water heat sink, or where the heating is extreme, 
a metal such as lithium can be used. Lithium's melting point is about 
3 6 0 0 F . and has a heat flux capacity of at least 1200 BTL/t t 2 -sec . 


69 





wback 


5 F as 
de an 



es a 
s e o u s 
id is 
f ue 1 . 
1 a 1 i er 
to t he 


ed for 
r a 1 i y , 

> . The 
tnd the 


70 




leading edge narrow the possibilities of systems capable of keeping it 
cool. 

The onlv type which cannot be used is transpiration. The high 
pressures at the leading edge and stagnation area preclude its use. 
While using some cooling system which ultimately uses the fuel to cool 
the aircraft appears attractive at this point, the fuel has a limited 
heat capacity which will quickly run low if used to cool the leading 
edge. Therefore the choice for leading edge cooling is a lithium spray 
cooling system which is relatively light( about 4QQlbs.) and self- 
contained. This system will provide cooling for 6" from the leading 
edge. Other systems are used elsewhere. 

The rest of the forward surfaces will need cooling, although not 
nearly as much as the leading edge. A transpiration system is desired 
here, but at this point it is unknown the effects it may have on the 
inlets and how it may be bled off before them. If transpiration is not 
possible, a direct convection system will be used to cool the forward 
surfaces. As the fluid is heated, it is transferred to the rear of the 
aircraft where it can be used for transpiration cooling on the exhaust 
expansion surfaces. If this "free" transpiration coolant is not 
available (i.e. transpiration is used on the forward surfaces), a metal- 
matrix surface will be used on the exhaust surfaces. 

Another critical area of heating is the propulsion area, including 
inlets. The large temperatures involved require an active cooling 
system. At this point a direct convection system looks appealing. A 
system of cooling plates will provide the necessary cooling via the 
hydrogen fuel as it is pumped to the engines. 
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The avionics bay, containing the communication and control 
electronics may need cooling to keep it below 10Q°F which may be 
required by these devices. If the fuel tank placement doesn’t provide 
the necessary cooling, an indirect convection system will be used. This 
system will use a water-glycol mixture to cool the avionics bay. A heat 
exchanger will be used with the fuel to utilize more of hydrogens heat 
capac i ty . 

The only areas not mentioned so far are the control surfaces. 
These surfaces can be exposed to extreme temperatures close to that of 
the forward surfaces. A metal-matrix material will be used on all 
control surfaces. 

A hypersonic aircraft will prove to be a flying cooling system as 
discussed above. The economical use of these technologies may well 
prove to be the decisive factor in the feasibility of a hypersonic 
ai rcraf t . 
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Chapter 12 


SUPPORT SYSTEMS 
(Jin Rogers ) 


Cont ro 1 

The control of a vehicle which flies at ten times the speed of 
sound has some interesting problems. At these speeds the temperature of 
the aircraft can adversely effect the control systems. To alleviate 
some of this problem, a fly-by-light system will be used. The system 
will be controlled by a dua 1 - redundant , f au 1 1 - to 1 eran t , centralized 
system located in the avionics bay. By using this centralized, 
integrated op t o - e 1 ec t ron i c system, the cooling of many areas may be 
eased . 

At this point, the effects of these temper atures on the hydraulic 
system is unknown. The hydraulic system is required to actuate the 
control surfaces, inlets/exits, and landing gear. This system, not 
infallible, requires a backup system. Under consideration is a 
duplicate system or a backup electrical system. 

The initial propulsion device will provide electrical power to the 
aircraft via generators. When the scramjet. is being tested, the 
generators will not be producing any power. The primary propulsion 
system will not be started until the aircraft has reduced its speed 
considerably in order to conserve fuel. Possibilities for this system 
range from a hydrazene power generator to a secondary turbine using a 
ram input flow, although there may be a possibility of using the ramjet 
turbine by diverting some of the airflow through it. The length of time 
an auxiliary power system will be operating will be anywhere from 4-12 
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minutes, depending on the final deceleration decided upon for the flight 
plan. 

Communi cat ion 

An unmanned aircraft can either have a programmed flight plan or it 
can be radio-controlled. This aircraft will have both. An aircraft 
flying at hypersonic speeds will require a nearly real-time interface 
with minimum delay between command and reaction. This requires the 
telemetry data to be received often and quickly in order to maintain a 
tight control on it. 

This requirement of quick communication requires a large bandwidth 
at high frequencies. As an arbitrary choice due to the uplink/downlink 
duality of the frequency range and the low signal absorbsion by the 
atmosphere, a X-band frequency will be used. Communicating with this 
aircraft will also require a number of stations. 

The aircraft will fly over 1S0Q miles in about 15 minutes. A 
proposed flight plan would have the aircraft launched from the carrier 
plane flying out of McChord A.F.B. outside of Seattle, Wa . 

The path then runs parallel to the west coast and will land at Edwards 
A.F.B. east of Los Angeles, Ca . A system of up 1 i nk/ down 1 i nk stations 
distributed along the coast will provide uninterrupted communication 
to/from the aircraft. 

The critical stages of the flight, as with any other, is the 
beginning and the end. The carrier/test plane combination will perform 
as an integrated system from takeoff through separation. This should 
provide a higher level of safety during this phase. During landing, a 


video camera is currently being considered to provide the controller a 
better feeling of the aircraft as it lands. 

This will all be for naught if something goes wrong in the 
communication chain, even a momentary lapse. By having a programmed 
flight plan in memory prior to separation, the test can be conducted 
without external control. To land, a backup chase plane will be able to 
control the test plane via different channels allowing the "backseat" to 
land it if all else fails. This method was developed for the Hi MAT 
research aircraft and seems re 

The control and communi- 
cation subsystems will become 
a very integrated, inter- 
locked system. Figure 12.1 
illustrates this system as it 
is currently envisioned. An 
interesting item is what type 
of antenna can be used at 
such high temperatures. 
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Figure 12.1: Control Diagram 
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Chapter 13 


COST ANALYSIS 
(Gary Hufford) 

The total cost of the research vehicle has been estimated to be 
3 . 7077 billion dollars. This figure was based upon statistical methods 
presented in Chapter 24 of "Fundamentals of Aircraft Design" by Leland 
Nicolai. This statistical method was based upon the top speed of the 
aircraft, the ampere weight of the aircraft, and the number of aircraft 
produced . 

The following numbers were used in the cost analysis 
computation. The ampere weight was estimated to be 15000 pounds. This 
weight is the empty weight minus weight due to engines, landing gear, 
avionics and cooling system. The top speed used was 5943 knots. This 
is the speed at Mach 10 and 100,000 feet. The number produced is one. 

A breakdown of the 
cost in eight different 
areas is given in Figure 
1 . The major areas of 
cost is in the number of 
engineering hours requ- 
ired and the developmen- 
tal support of the vehicle. These are expected to be the major costs of 
a research vehicle of this type. The carrier aircraft was not figured 
into the costs of the research vehicle. The cost of the carrier 
aircraft is important and we have not overlooked designing the aircraft 
to minimize the cost of the carrier aircraft. Our small design will 


Engineering Hours 

1 ,880 ,000 ,000 

Developement Support 

1 , 175,000,000 

Flight Test Operations 

30,500, 000 

Too 1 ing 

403 ,800 ,000 

Manufacturing Labor 

93 , 930 , 000 

Quality Control 

122 ,000 , 000 

Material and Equipment 

1,030, 000 

Engine Costs 

1 , 440 . 000 

Total Cost for 1 Aircraft 

= 3.7077 Billion 


Figure 13.1: Cost Analysis 
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help minimize the cost of the carrier aircraft, possibly even eliminat- 
ing the need for a new design if our design can possibly be launched 
from under the wing of a B-52. 


CONCLUSION 

(Gary Hufford} 

We have shown that the vehicle described in this report can test 
scramjet engines up to a maximum Mach number of 9.9 at a cost approxi- 
mately 3.7 billion dollars. 

We feel that this design has three major advantages. First, it is 
a very small design. A small design is significantly cheaper when the 
carrier aircraft is figured into the cost analysis. A small design also 
minimizes the sonic boom problem present with a supersonic aircraft. 
Second, the aircraft can be used in an efficient manner to test scramjet 
engines at low Mach numbers ( 6 to 8 ) as well as higher Mach numbers { 
9 to 10 ). The low Mach number tests can be conducted with the base 
design and the high Mach number tests can simply add the rocket 
attachments to attain higher Mach numbers. The aircraft is not designed 
for a specific Mach number, rather it is a general design for a 
hypersonic airplane. Therefore, it does not require a specific Mach 
number for optimization. The aircraft is versatile in the selection of 
Mach numbers that the scramjets can be tested. Third, the aircraft has 
a separate engine system ( turbof anramje t ) that can be used if the 
scramjet tests fail. The t u r bo f anramje t uses the same fuel as the 
scramjets. Thus, if the scramjet fails, the fuel meant for the scramjet 
can be diverted to the t urbo f anr am j e t to land the vehicle safely. 
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Figure A. 4: Configuration 4 
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Figure A. 5: Configuration 5 
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GE Turbof cxnram jet 




Figure A. 6: Engine location 
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GE Turbof anram jet 



37 ' 



Figure A. 7: Reference values 
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SCALE; 








Ing Loading 44.7 Ib/ft^ 



Figure A. 9: Final configuration 
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GE Turbof anram jet 



Appendix B 

Weights Analysis Program 


40000.00000 ESTIMATED GROSS WEIGHT EGW 

2905.0 DELTA WING WT COEF WINGC 

3.0 ULTIMATE LOAD FACTOR ULF 

16.0 STRUC. SPAN (ACROSS MID CHORD) SPAN 

205.0 WING PLANFORM AREA SWING 

0.65 WING ROOT THICKNESS ROOT 

0.621 DELTA WING WT EXP WINGE 

3.9 VERT TAIL WT COEF VTC 

400.0 PLANFORM AREA SVERT 

1.00 WT EXPONENT VTE 

1000.0 BODY WETTED AREA SWET 

60.0 BODY LENGTH LBOD 

1500.0 MAX. DYNAMIC PRESSURE (Qmax) QMAX 

5.87 BODY HEIGHT HBOD 

1.0 BODY WT. EXP. BODE 

1.275 BODY TEMP WT MOD FACT (500 F) BTMOD 

33000.0 MAX TURBOFANRAMJET NET THRUST TTURBO 

5000.0 MAX SCRAMJET NET THRUST TSCRAM 

0.00625 THRUST STRUCT. WT COEFF TSC 

6 9.0 WT INTERCEPT TSI 

0 . 4 SECONDARY BODY WT COEF SBODC 

2000.0 THERM PROTECTION SURFACE AREA STPS 

0.7 BODY INSUL. UNIT WEIGHT COEF BINSC 

1.4 BODY COVER PANELS UNIT WT COEF BCOVC 

0.0025 LAUNCH GEAR WT. COEFF LGC 

0.00916 AIR LAUNCHED LNDG GEAR WT COEF LNDGC 
1.124 WT EXP LNDGE 

5400.0 INDIVIDUAL TURBORAMJET WEIGHT TURBO 

1580.0 INDIVIDUAL SCRAMJET WEIGHT SCRAM 

1.0 # OF TURBORAMJETS XTURBO 

3.0 # OF SCRAMJETS XSCRAM 

0.004 ENGINE MOUNT WT. COEFF. ENGMTC 

0.53 LH2 INTEGRAL TANK WT COEF TANKC 

1932.0 FLIGHT TIME (SEC) TFLIGHT 

400.0 FUEL TANK RADIATING TEMP (F) TRAD 

4.345 DUCT (INTERNAL) WT INTERCEPT DII 

4.0 LIP TO ENGINE FACE LENGTH (FT) LINLET 

2.0 NUMBER OF INLETS XINLET 

5.8 INLET CAPTURE AREA ( FT~ 2 ) SCAPT 

118.0 INLET PRESSURE (PSIA) PINLET 

9.0 DESIGN MACH NUMBER DMACH 

117.35 INLET RAMP WT COEF RAMPC 

8.0 RAMP LENGTH (FT) LRAMP 

0.294 INLET RAMP WT EXP. RAMPE 

35.0 INLET SPIKE WT COEF SPIKEC 

0.323 AERO ORIENT. CONTROL WT COEF CONTC 

0.903 EXP CONTE 

37^0 PLANE GEOMETRIC SPAN (FT) GSPAN 
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1.167 

4.64 

30.0 

0.361 

0.018 

0.006 


7100.0 

2000 . 0 
5 


3.1 

2.45 

3 . 1 

2.45 

1.3 

1.47 

1.3 

1.47 

4.08 

0.98 

CAPACITY 

7000 . 0 


ELECTRIC. SYS WT COEF . ELECC 

HYDRAULIC SYS WT COEF. HYDRC 

AVIONICS SYS WT COEF AVIONC 

EXP. AVI ONE 

TRAPPED FUEL WT. COEF TRFUC 

INFLIGHT FUEL LOSS COEF FLOSSC 

WEIGHT OF USEABLE FUEL WUFUEL 


ACTIVE COOLING 
NUMBER OF FUEL 

368.0 

368.0 

170.0 

170.0 

456.0 

ROCKET BOOSTER 


WT . ALLOWANCE ACOOL 

TANKS NTANKS 

TANKA ( I ) , TANKB ( I ) , VTANK ( I ) 

TANKA= FUEL TANK MAJOR DIA 
TANKB = MINOR DIA 

VTANK = FUEL TANK VOLUME 

WEIGHT ROCKETS 
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PROGRAM WEIGHT 


★ 

* 

* 

★ 

* 

* 

* 

* 

* 

Hr 


This program is based on the Weights Analysis of Advanced 
Transportation Systems (WAATS) written by C.R. Glatt. NASA 
Report CR-2420. 


The program takes formatted input from a data file that has a 
description of each variable following the 

numerical value. This allows input files to be edited easily 
on a word processor by changing the appropriate numeric values. 


REAL LBODY , LGC , LNDGC , LNDGE , LRAMP , L INLET 


Hr 

*************** data input 


Hr 


OPEN ( UNIT=1, FI LE=' WE IGHT.DAT' , STATUS= ' OLD ' ) 


READ (1,1) 

EGW 

READ (1,1) 

WINGC 

READ (1,1) 

ULF 

READ (1,1) 

SSPAN 

READ (1,1) 

SWING 

READ (1,1) 

ROOT 

READ (1,1) 

WINGE 

READ (1,1) 

VTC 

READ (1,1) 

SVERT 

READ (1,1) 

VTE 

READ (1,1) 

SWET 

READ (1,1) 

LBODY 

READ (1,1) 

QMAX 

READ (1,1) 

HBODY 

READ (1,1) 

BODE 

READ (1,1) 

ETMOD 

READ (1,1) 

TTURBO 

READ (1,1) 

TSCRAM 

READ (1,1) 

TSC 

READ ( 1 , 1) 

TSI 

READ (1,1) 

SBODC 

READ ( 1,1) 

FTPS 

READ (1,1) 

BINSC 

READ (1,1) 

P.COVC 

READ ( 1,1) 

I.GC 

READ ( 1,1) 

LNDGC 

READ ( 1,1) 

LNDGE 

READ ( 1,1) 

TORBO 

READ ( 1,1) 

SCRAM 

READ ( 1,1) 

XTURBO 

READ ( 1,1) 

XoCRAM 

READ ( 1,1) 

ENGMTC 

READ ( 1,1) 

TANKC 

READ ( 1,1) 

TFLIGHT 
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READ (1,1) 

TRAD 

READ (1,1) 

FTSWET 

READ (1,1) 

DII 

READ (1,1) 

L INLET 

READ (1,1) 

XINLET 

READ (1,1) 

SCAPT 

READ (1,1) 

PINLET 

READ (1,1) 

DMACH 

READ (1,1) 

RAMPC 

READ (1,1) 

LRAMP 

READ (1,1) 

RAMPE 

READ (1,1) 

SPIKEC 

READ (1,1) 

CONTC 

READ (1,1) 

CONTE 

READ (1,1) 

GSPAN 

READ (1,1) 

ELECC 

READ (1,1) 

HYDRC 

READ (1,1) 

AVIONC 

READ (1,1) 

AVI ONE 

READ (1,1) 

TRFUC 

READ (1,1) 

FLOSS 

READ (1,1) 

WUFUEL 

READ (1,1) 

ACOOL 

CLOSE (UNIT=1) 


FORMAT (FI 2. 5) 



* GENERAL CALCULATIONS 
WTO= EGW 

* Max total thrust 

TTOT= XTURB.O*TTURBO + XSCRAM*TSCRAM 

* Temperature correction factor 

TEMPFCT= 0.203 * DMACH + 0.4 

* Gross fuel at takeoff 

GRFUEL = WUFUEL+ FLOSS *WUFUEL 

* Fuel tank volume to hold gross fuel (LH2) 

VTANK= GRFUEL/ 4.389 

* Landing weight 

10 WLAND= WTO - GRFUEL 

* 

*************** weight calculations 
* 

* weight of aerodynamic surfaces 

WWING= WINGC* ( WTO*ULF*SS PAN* SWING/ROOT ) **WINGE* 1 . 0E - 6 
WVERT= VTC*SVERT**VTE 
WAERO= WWING + WVERT 

* 

* weight of body structure 

* 

A= 0.341 * BTMOD 

B= (LBODY*ULF/HBODY)**0. 15 * QMAX**0.16 * SWET**1.05 


* 

* 

* 


* 

it 

* 


it 

it 

it 


it 

it 

it 


it 

★ 

* 


* 

★ 


it 


WBASIC= A*B**BODE 
WTHRST= TSC * TTOT + TSI 
WSECST= SBODC * SWET 
WSTRUCT= WBASIC + WTHRST + WSECST 

thermal protection 

WINSUL= BINSC * STPS 
WCOVER= BCOVC * STPS 
WTPS= WINSUL + WCOVER + ACOOL 

launch and landing gear 


WLAUNCH= LGC * WTO 

WLG= LNDGC * WLAND * * LNDGE 

WGEAR= WLAUNCH + WLG 

engine system weight 


WTURBO= XTURBO * TURBO 

WSCRAM= X SCRAM * SCRAM 

WENGMT= ENGMTC * TTOT 

WENGINES= WTURBO + WSCRAM + WENGMT 


inlet system weight 

B=1 . 733* ( L INLET* XI NLET) ** . 5* ( SCAPT/XINLET) ** . 3334*PINLET** . 6667 
wTnnrT= dii * B 

WRAMP= RAMPC* ( LRAMP*XINLET* ( SCAPT/XINLET) ** . 5*TEMPFCT ) **RAMPE 
WSPIKE= SPIKEC * XTURBO 
WINLET= WIDUCT + WRAMP + WSPIKE 

fuel system 


WFUNCT= TANKC * VTANK 

TCORR= 0.7 + TFLIGHT*8 . 0E-5 
A= TCORR* ( 0 . 0007 * TRAD + 0.24) 

WINSFT= A * FTSWET 

A= 0.00651 - 2.995E-8 * TTOT 
B= .1 + .175 * LBODY 
WFUSYS= A*TTOT + B*LBODY 
WCRYO= WFUNCT + WINSFT + WFUSYS 

controls and avionics 

A= WTO**0.6667*(LBODY + GSPAN)**0.25 
WORIENT= CONTC * A**CONTE 

WELECT= ELECC * (SQRT(WTO) * LBODY**0.25) 

C= ( LBODY+SSPAN ) **0 . 5 * 1.25 
WHYDR=HYDRC* ( { SWING+SVERT ) *QMAX/ 1000) **.334 + C 
WAVON= AVI ONC* WTO** AVI ONE 
WCONT= WORIENT + WELECT + WHYDR + WAVON 
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final summation 


WNEW= WAERO+WSTRUCT+WTPS+WGEAR+WENGINES+WINLET+WCRYO+WCONT 

&+GRFUEL 


iteration loop 

ERR= ABS ( WTO-WNEW) 

IF ( ERR .LT. 0.01 ) GOTO 20 
WTO= WTO* ( WNEW/ WTO ) * * 1 . 1 
GOTO 10 


print output 


OPEN ( UNIT= 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE ( 2 , * ) 
WRITE ( 2 , * ) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE { 2 , * ) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2 , * ) 
WRITE ( 2,*) 
WRITE ( 2 , * ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2,* ) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE ( 2 , * ) 


2, FILE=' WEIGHT. OUT' , STATUS='NEW' ) 

• FINAL WEIGHT ESTIMATES' 


'Gross takeoff weight ',WTO 

'Empty weight ' > WLAND 

'Gross fuel weight ' , GRFUEL 

'Useable fuel weight ' ,WUFUEL 

' Aerodynamic Surfaces ',WAERO 

'Wing structure ',WWING 

'Vertical stabilizer ' ,WVERT 

' Body Structure ' ,WSTRUCT 

'Main body ' ,WBASIC 

'Secondary body structure ',WSECST 

'Thrust structure ' , WTHRST 

' Environmental Protection ' ,WTPS 

'Structure insulation ' ,WINSUL 

'Body cover panels ' ,WCOVER 

'Active cooling allowance ',ACOOL 

' Launch & Landing Gear ' , WGEAR 

'Launch gear ' r WLAUNCH 

'Landing gear ',WLG 

' Engine System ' ,WENGINES 

' Turbof anram j et ' ,WTURBO 

'Scramjet ' ,WSCRAM 

'Engine mounts ' ,WENGMT 

' Inlet System ' , WINLET 

'Internal ducts ' ,WIDUCT 

'Variable ramps ' ,WRAMP 

' SDi ke ' , WSPIKE 


' Cryogenic Fuel System 
'Integral fuel tank 


' , WCRYO 
' , WFUNCT 




WRITE (2,*) 'Fuel tank insulation 
WRITE ( 2 , * ) 'Fuel distribution system 
WRITE (2,*) 

WRITE(2,*) ' Control and avionics system 

WRITE(2,*j 'Aerodynamic orientation controls 

WRITE ( 2 , * ) 'Electrical system 

WRITE (2,*) 'Hydraulic system 

WRITE ( 2 , * ) 'Avionics 

CLOSE ( UNIT= 2 ) 

STOP 

END 


' , WINSFT 
' , WFUSYS 

' , WCONT 
' , WORIENT 
' , WELECT 
' , WHYDR 
' , WAVON 
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QUO boo 


Appendix C.l 
Subsonic program 


IMPLICIT REAL (A-Z) 

DIMENSION T ( 205 ) , RHO(205), P(205), 0(205,105) 

INTEGER I , J , MSUB , ALTSUB , SELE 

OPEN(1 ,FILE=’A515.DAT’ ,STATUS=’NEW’ ) 

OPEN ( 4 , FILE= ’ GECM1 .DAT’ ,STATUS=’OLD’ ) 

OPEN (12, FILE= ’ NICOL . DAT ’ ,STATUS=’OLD’ ) 

OPEN ( 2 , FILE= ’ CL . DAT ’ , STATUS* ’ NEW’ ) 

OPEN ( 3 , FILE= ’ CD. DAT ’ , STATUS* ’NEW'-) 

OPEN ( 7 , FILE= ’ MACH . DAT ’ ,STATUS='NEW’ ) 

OPEN ( 8 , FILE= ’ ALT . DAT ’ , STATUS* ’ NEW’ ) 

OPEN (9, FILE* ’LD. DAT’ .STATUS* ’NEW’ ) 

OPEN( 10, FILE*’ DRAG. DAT’ , STATUS* ’ NEW’ ) 

OPEN (11, FILE* ’LIFT. OAT’ , STATUS* ' NEW’ ) 

OPEN (13, FILE* ’ ALPHA . DAT ’ , STATUS* ’NEW’ ) 

R = 1716. 

G0 = 32.2 
PI = 3.14157 
GAM = 1 .4 


************************ 

Atmosphere from 0 -> 36000 ft 


DO 6 1=1 ,73 
PI = 2116.2 
T1 = 518.69 
A = -.00356624 
H = 500* (I- 1 ) 

T< I ) = T1+A-H 

P(I) = P1*EXP(-G0/(R*T(I) )*H) 
RHO(I) = P(I)/(T(I)*R) 

6 CONTINUE 


Atmosphere from 36500 -> 82000 ft 


DO 7 1=74,165 
T 1 = 389.99 
PI = 464.86 
HI = 36500 
H = 500* (1-1 ) 
T( I ) = T1 


ooo oooo ooo 


P(I) = PI *EXP( -G0/ (R*T1 )* (H-H1 ) ) 
RHO(I) = P(I)/(T1*R) 

7 CONTINUE 


Atmosphere from 82500 -> 100000 ft 


DO 10 1=166,201 
HI = 82500 
PI = 55.837 
T1 = 390.24 
A = .001646 
H = 500* (1-1 ) 

T(I) = (H-H1 )*A+T1 

P(I) = PI *EXP( - ( (G0/ (R*T(I) ) )* (H-H1 ) ) ) 
RHO(I) = P(I)/ (T(I)*R) 

10 CONTINUE 


Dynamic pressures at various altitudes and mach numbers 
Q(alt,M) 


DO 15 ALTSUB= 1,201 
DO 20 MSUB=1 , 101 
M = 0.1* (MSUB- 1 ) 

Q(ALTSUB.MSUB) = P(ALTSUB)*(GAM/2)*M*M 
20 CONTINUE 
15 CONTINUE 


Airflow Data 


1000 WRITE(6,*) (’ Starting Mach number (.1 min):’) 
READ(5, * ) SM 

WRITE ( 6 , * ) (’ Ending Mach number (.9 max):’) 
READ<5,*) EM 

WRITE ( 6 , * ) (’ Step size (.1 multiples):’) 
READ(5 , * ) STM 


WRITE ( 6 , * ) (’ Starting altitude (0 ft min):’) 
READ(5, * ) SALT 

WRITE ( 6 , * ) (’ Ending altitude (100000 ft max):’) 
READ(5, * ) EALT 

WRITE ( 6 , * ) (’ Step size (500 ft multiples):’) 
READ(5,* ) STALT 


WRITE (6, * ) (’ Starting angle of attack:’) 
READ(5 , * ) SALP 

WRITE (6, * ) (’ Ending angle of attack:’) 


o o o o o 


READ(5,*) EALP 

WRITE ( 6 , * ) <’ Step size:’) 

READ(5,* ) STALP 

WRITE ( 6 , * ) (’ Max drag:') 
READ(5, * ) MAXD 

WRITE ( 6 f * ) (’ Max lift:’) 
READ(5,*) MAXL 

WRITE ( 6 , * ) (’ Min lift:’) 
READ(5, * ) MINL 


c ***************** 
C Geometry Data 
c ****»************ 


WRITE ( 6 , * ) (’ Span (ft):’) 

READ(4, * ) SPAN 

WRITE ( 6 f * ) (’ Wing ref. area SE (ft*2):’) 
READ(4, * ) SE 

WRITE ( 6 f * ) ( ’ Wetted wing area (ft*2):’) 

READ ( 4 , * ) SWETW 

WRITE ( 6 , * ) (’ Wetted body area (ft*2):’> 

READ(4 , * ) SWETB 

WRITE ( 6 , * ) (’ Max wing length (ft):’) 

READ(4 , * ) WL 

WRITE ( 6 , * ) (’ Max body length (ft):') 

READ(4 , * ) BL 

WRITE ( 6 , * ) (’ Sweep angle (deg):’) 

READ(4 , * ) SAP 

WRITE ( 6 , * ) (’ Max wing thickness to chord:’) 
READ(4,*) TC 

WRITE ( 6 , * ) (’ Leading edge radius (% chord):’) 
READ(4,* ) RLE 

WRITE(6,*) (’ Body cross-section area (ft“2):') 
READ(4,*) SB 

WRITE ( 6 , * ) (’ Body width (ft):’) 

READ(4,*) BD 


9WPR = SWP* (PI/ 180) 
AR = ( SPAN*SPAN ) / SE 


Subsonic Inputs 

If changing airflow properties only, 

the proceeding variables do not need to be changed. 


DBD = BD/SPAN 
CSWP = COS(SAPR) 


bob b o o o 


WRITE (6,100) DBD 

100 FORMAT ( 1 X , ’ Get F for d/b = ’,F7.5,’ from FIG. 11.3: ) 

READ( 12, * ) F 
WRITE (6,102) AR 

102 FORMATdX,’ Get Cl for AR « ’ ,F12.5,’ from FIG. 2.6: ) 

READ( 12, * ) Cl 

WRITE ( 6 , 1 06 ) AR , _ 

106 FORMATdX,’ Get e' for AR = ’,F12.5,’ from FIG. 11.5: ) 

READ(12, * ) EPR 

WRITE (6,108) RLE ^ t4 

108 FORMATdX,’ Get K’ ’ for RLE = ’,F12.5,’ from FIG. 11.6: ) 

READ(12, * ) KDPR 

DO 55 MACH=SM,EM,S7M 

DO 50 ALT=SALT , EALT , STALT 

DO 60 ALPHA D=SALP , EALP , STALP 
ALPHA = ALPHAD* ( PI / 1 80 ) 

ALTSUB » ALT/ 500+1 
MSUB = MACH/. 1+1 
MAC = .5*WL 

MU = . 317E- 10* (T( ALTSUB )**1 .5* (734.7/ (T(ALTSUB)+216) ) ) 
B = SORT d -MACH* *2) 


Subsonic Lift 


CLAW/1 = 1 + (TAN(SWJPR) ) * *2/ (B*B) 
CLAW/2 = SORT ( 4+AR*AR*B*B*CLAW/1 ) 
CLAW/ = (2*PI*AR) / <2+CLAW/2) 

CLAWJB = F*CLAW 

CL = CLAW8*ALPHA+C1 ‘ALPHA* *2 


Subsonic Drag 
CD0w 


IF (MACH. LT. 0.25) THEN 
RR= . 7407*CSW/P+ . 43072 
ENDIF 

IF ( (MACH. GT. 0.25) .AND. (MACH. LT. 0.6) ) THEN 
RR= . 7407*CSWP+ . 2086*MACH+ . 37857 
ENDIF 

IF ( (MACH. GT. 0.6) .AND. (MACH. LT. 0.8) ) THEN 
RRi . 7407*CSWR+ . 635*MACH+ . 1 2272 
ENDIF 

IF ( (MACH. GT. 0.8) .AND. (MACH. LT. 1 ) ) THEN 
RR=.7407*CSWR+1 . 18*MACH- .31328 
ENDIF 


REW/ = ( RHO( ALTSUB ) *MAC*MACH* SORT (GAM*R*T (ALTSUB) ) )/NU 
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CFP = .455/ (LOG10(REW) )**2.58 

CD0W = CFP*( 1+1 . 2*TC+100*TC**4)*RR*(SWETW/SE) 


CD0B 


DE * SORT (SB/. 7854) 

REB = (RHO(ALTSUB)*BL*MACH*SQRT(GAM*R*T(ALTSUB) ) )/MJ 

CF = .455/ (LOG10( REB) ) **2.58 

CDFB = CF* (1+60/ (BL/DE)**3+.0025* (BL/DE) )* (9AIETB/SE) 
CD0B = CDFB 
CD0 = CD0W*CD0B 


K' 


E = EPR* ( 1 - ( BD/ SPAN ) * *2 ) 
KPR = 1 / (PI*AR*E) 


K’ ’ 


K = KPR+KDPR 
CD = CD0+K*CL**2 


DRAG = CD*Q(ALTSUB, MSUB)*SE 
LIFT = CL*Q(ALTSUB, MSUB)*SE 
LD = CL/CD 

IF ( (LIFT.GT.MINL) .AND. (LIFT. LT.MAXL) .AND. 

* ( DRAG . LT . MAXD) ) THEN 

WRITE (1,1 10) MACH, ALT, ALPHAD 

110 FORMAT (IX,’ MACH: \F7.3,’ ALT: ',F12.4,' ALPHA: *,F7.3) 
WRITE(1,112> LIFT, DRAG 

112 FORMAT (IX,’ LIFT: *,F14.7,’ DRAG: ’.F14.7) 

WRITE(1,114) CL, CD 

114 FORMAT (IX,’ CL: ',F8.5,’ CD: ’,F8.5) 

WRITE ( 1,116) Q ( ALTSUB , MSUB ) , LD 
116 FORMAT (IX,’ DYNAMIC PRESSURE= *,F12.5,' L/D: *,F7.4) 

write(i,*) r ’) 

WRITE (2, * ) CL 
WRITE (3, * ) CD 
WRITE(7,*) MACH 
WRITE ( 8 , * ) ALT 
WRITE ( 9 , * ) LD 
WRITE ( 10, * ) DRAG 
WRITEdl,*) LIFT 


1 Q 3 


WRITE < 13,*) ALPHAD 
ENDIF 

60 CONTINUE 

50 CONTINUE 
55 CONTINUE 

1008 WRITE ( 6 , * ) (’ Enter Selection: ’ ) 

WRITE ( 6 , * ) (’ 1. Change geometry only’ ) 

WRITE (6, * ) (’ 2. Change airflow properties only’) 

IF (SELE.EQ.1) THEN 
GOTO 1000 
ENDIF 

STOP 

END 
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Appendix C.2 
Supersonic program 


IMPLICIT REAL (A-Z) 

DIMENSION T ( 205 ) , RHO(205), P(205), (3(205,105) 

INTEGER I , J , MSUB , ALTSUB , FLAG , SELE 

CHARACTER *8 FCL, FCD, FCLD, FM, FALT 

OPEN ( 1 , FILE= ’A515.DAT’, STATUS= » NEW’ ) 

OPEN ( 4 , FILE® ’ GEOM1 . DAT ' , STATUS= » OLD’ ) 

OPEN ( 2 , FILE® ’ CL . DAT ’ .STATUS® ’NEW’ ) 

OPEN ( 3 , FILE= ’ CD . DAT ’ , STATUS= ’ NEW* ) 

OPEN (7 , FILE= ’ MACH .DAT’ , STATUS® ' NEW’ ) 

OPEN ( 8 , FILE® ’ ALT . DAT ’ , STATUS® 'NEW' ) 

OPEN (9, FILE®’ LD. DAT’ .STATUS® ’NEW’ ) 

OPEN (10, FILE® ’ DRAG . DAT ’ .STATUS® ’NEW’ ) 

OPEN (11, FILE® ’ LIFT . DAT ’ , STATUS® ’ NEW’ ) 

OPEN (12, FILE® ’ ALPHA . DAT ’ , STATUS® ’ NEW’ ) 

R = 1716. 

G0 = 32.2 
PI = 3.14157 
GAM = 1 .4 


Atmosphere from 0 -> 36000 ft 


DO 6 1=1,73 
PI = 2116.2 
T1 = 518.69 
A = -.00356624 
H = 500* (1-1 ) 

T( I ) = T1+A*H 

P(I) = PI *EXP( -G0/ (R*T( I) ) *H) 
RHO(I) = P(I) / (T(I)*R) 

6 CONTINUE 


Atmosphere from 36500 -> 82000 ft 


DO 7 1=74,165 
T 1 = 389.99 
PI = 464.86 
HI = 36500 
H = 500* (1-1 ) 



T(I) = T1 

P(I) = P1*EXP(-G0/(R*T1 )*(H-H1 ) ) 
RHO(I) = P(I)/(T1*R> 

7 CONTINUE 


C Atmosphere from 82500 -> 100000 ft 
c ******************************** 

DO 10 1=166,201 
HI = 82500 
PI = 55.837 
T1 = 390.24 
A = .001646 
H = 500* (I-D 
T(I) = (H-H1 ) *A+T1 

P(I) = P1*EXP(-((G0/(R*T(I)))*(H-H1))) 
RHO(I) = P(I)/(T(I)*R) 

10 CONTINUE 


C Dynamic pressures at various altitudes and mach numbers 
C Q(a1t,M) 

c ***********************«*****«********** 

DO 15 ALTSUB=1 ,201 
DO 20 MSUB=1 ,101 
M = 0.1* (MSUB-1 ) 

Q ( ALTSUB , MSUB ) = P(ALTSUB)*(GAM/2)*M*M 
20 CONTINUE 
15 CONTINUE 


C Airflow Data 


1000 WRITE (6,*) (’ Starting Mach number (1.3 min):’) 
READ(5, * ) SM 

WRITE(6,*) (’ Ending Mach number (5 max):’) 
READ(5,*) EM 

WRITE ( 6 , * ) (’ Step size (.1 min):’) 

READ(5, * ) STM 


WRITE ( 6 , * ) (’ Starting altitude (0 ft min):’) 
READ (5 *) SALT 

WRITE(6,*) (’ Ending altitude (100000 ft max):’) 
READ(5, * ) EALT 

WRITE ( 6 , * ) (’ Step size (500 ft min):’) 

READ ( 5 , * ) STALT 


WRITE<6, * ) (’ Starting angle of attack:’) 
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READ(5, * ) SALP 

WRITE ( 6 , * ) (’ Ending angle of attack:’) 

READ(5,*) EALP 

WRITE ( 6 , * ) (’ Step size:’) 

READ(5, * ) STALP 

WRITE ( 6 , * ) (’ Max drag:’) 

READ<5,*) MAXD 

WRITE ( 6 , * ) (’ Max lift:’) 

READ(5, * ) MAXL 

WRITE ( 6 , * ) (’ Min lift:’) 

READ(5,*) MINL 


C Geometry Data 

c **********«***************************** 

WRITE ( 6 , * ) (’ Span (ft):') 

READ(4, * ) SPAN 

WRITE ( 6 , * ) (’ Wing ref. area SE (ft*2):’) 

READ(4, * ) SE 

WRITE ( 6 , * ) ( ’ Wetted wing area (ft*2) ’) 

READ(4 , * ) SWETW 

WRITE ( 6 , * ) (’ Wetted body area (ft"2):’) 

READ(4 , * ) SWETB 

WRITE ( 6 , * ) (' Max wing length (ft):’) 

READ(4, * ) WL 

WRITE ( 6 f * ) (’ Body length (ft):’) 

READ(4,*) BL 

WRITE ( 6 , * ) (’ Sweep angle (deg):') 

READ(4, * ) SAP 

WRITE ( 6 , * ) (’ Max wing thickness (% chord):') 
READ( 4 , * ) TC 

WRITE ( 6 , * ) (’ Leading edge radius (% chord):’) 
READ(4,* ) RLE 

WRITE ( 6 , * ) (’ Body cross-section area (ft*2):’> 
READ(4,*) SB 

WRITE ( 6 , * ) (’ Body width (ft):’) 

READ(4,*) BD 


* Supersonic 


DO 55 MACH=SM , EM , STM 
SWPR = SWP* (PI/ 1 80) 
B = SORT (MACH* *2- 1 ) 
AR = (SPAN* SPAN )/SE 
MAC = .5*WL 
ATS = AR*TAN(SWPR) 
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TSW = TAN(SWPR)/B 
SWET = SAETB+SAETW 


IF (TSW.LE.1) THEN 

WRITE ( 6 , 208 ) ATS, TSW 

208 FORMAT (IX,* When AR Tande sweep) =’,F12.5,’ and 

* tande)/B=’ ,F12.7, * get Cna from FIG. 11.2:’) 

ELSEIF (TSW.GT. 1 ) THEN 
TSW = 1 /TSW 
WRITE ( 6 , 210) ATS, TSW 

210 FORMAT( IX, ’ When AR Tande sweep) =',F12.5,’ and 

* B/TAN(LE)=’ ,F12.7, * get Cna from FIG. 11.2:’) 

ENDIF 

READ(5, * ) BCLA 

DO 50 ALT =SALT , EALT , ST ALT 

DO 60 ALPHAD=SALP , EALP , STALP 
ALTSUB = ALT/ 500+1 
MSUB = MACH/ .1+1 
ALPHA = ALPHAD* ( PI / 1 80 ) 

MU = .317E-10*(T(ALTSUB)**1 .5* (734.7/ (T(ALTSUB)+216) ) ) 
CLAWB = BCLA/B 
CL = CLAWB* ALPHA 


C *********************« — ** — 

C SUPERSONIC DRAG 

C CD0 OF THE WING (SUPERSONIC) 

C SHARP NOSED SECTIONS IE. LEADING EDGE RADIUS 
c LESS THAN .1% T/C 

c ***. ************************* 


IF (BCS.GE.1) THEN 

CDW = ( 16/3) /B*TC**2 
ELSEIF (BCS.LT.1) THEN 

COW = ( 16/3)* ( 1 /TAN(SWPR) )*TC**2 
ENDIF 

REB = ( RHO( ALTSUB ) *BL*MACH*SQRT (GAM*R*T (ALTSUB ) ) ) /MJ 
CF = .455/ (LOG10(REB) )**2.58 
CFCCFI = 1 / ( 1+. 144*MACH**2)**.65 
CDFB = CFCCFI*CF*SWETB/SE 

REW = ( RHO ( ALTSUB ) *WL*MACH* SORT (GAM*R*T (ALTSUB) ) ) /MJ 
CDFW = CFCCFI*CF*3A£TW/SE 
CD0 = CDFB+CDFW+CDW 


c ****. *****************«»---"'-- 

C DRAG DUE TO LIFT; SUPERSONIC 

c *. ****.*********************•* 


K = 1 /CLAWB 

CD = CD0 + K*CL**2 



110 

112 

114 

116 


60 

50 

55 


DRAG = CD*Q(ALTSUB, MSUB)*SE 
LIFT = CL*Q(ALTSUB, MSUB)*SE 
LD = CL/CD 

IF ( (LIFT. GT.MINL) .AND. (LIFT. LT.MAXL) .AND. 
(DRAG.LT.MAXD) ) THEN 

WRITEd,*) (’ ’> 

WRITE (1,110) MACH, ALT, ALPHAD 
FORMATdX,’ MACH: ’,F7.3,' ALT: ',F12.4, 

*’ ALPHA: ’ , F7 . 3 ) 

WRITE (1,112) LIFT, DRAG 

FORMAT (1 X , ’ LIFT: \F14.7,’ DRAG: *,F14.7) 
WRITE (1,114) CL, CD 

FORMAT ( 1 X , ’ CL: ',F8.5,’ CD: ’,F8.5) 
WRITE (1,116) Q(ALTSUB,MSUB) , LD 
FORMAT (1 X , ’ Q= *,F 12.6,’ L/D= »,F10.6) 

WRITE(2,*) CL 
WRITE ( 3 , * ) CD 
WRITE ( 7 , * ) MACH 
WRITE ( 8 , * ) ALT 
WRITE(9,*) LD 
WRITE(10,*) DRAG 
WRITEdl,*) LIFT 
WRITE (12,*) ALPHAD 
ENDIF 
FLAG=1 

CONTINUE 

CONTINUE 

CONTINUE 

WRITE ( 6 , * ) (’ Enter Selection:') 

WRITE(6,*) (' 1. Change airflow’ ) 

WRITE ( 6 , * ) (’ 2. Quit’) 

READ(5,*) SELE 
IF (SELE.EQ. 1 ) THEN 
GOTO 1000 
ENDIF 


STOP 

END 
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Appendix C.3 
Hypersonic program 


IMPLICIT REAL (A-Z) 

DIMENSION T ( 205 ) , RHD(205), P(205), 0(205,105) 

INTEGER I, J,MSUB,ALTSUB, SELE, GUIDON 
CHARACTER* 8 QUEST 

OPEN( 1 ,FILE=’A515.DAT’ , STATUS* * NEW' ) 

OPEN ( 4 , FILE= ’ GECM1 . DAT ' .STATUS* 'OLD' ) 

OPEN ( 2 , FILE= ’ CL . DAT ’ .STATUS* 'NEW' ) 

OPEN ( 3 , FILE= ’ CD . DAT ’ , STATUS= ’ NEW' ) 

OPEN ( 7 , FILE= ’ MACH .DAT’ , STATUS= ' NEW' ) 

OPEN (8, FILE*’ ALT . DAT ’ , STATUS= ' NEW' ) 

OPEN ( 9 , FILE= ’ LD . DAT ’ , STATUS* ’ NEW' ) 

OPEN ( 10, FILE* 'DRAG. DAT’ ,STATUS*’NEW’ ) 

OPEN (11, FILE* ’ LIFT . DAT ’ , STATUS* ’ NEW’ ) 

OPEN ( 1 2 , FILE* ’ ALPHA . DAT ’ , STATUS* ’ NEW’ ) 

R = 1716. 

G0 = 32.2 
PI = 3.14157 
GAM = 1 .4 


************************ 

Atmosphere from 0 -> 36000 ft 


DO 6 1=1,73 
PI = 2116.2 
T1 = 518.69 
A = -.00356624 
H = 500* (1-1 ) 

T( I ) = T1+A*H 

P(I) = PI *EXP( -G0/ (R*T ( I) ) *H) 
RHO(I) = P(I)/(T(I)*R) 

6 CONTINUE 


Atmosphere from 36500 -> 82000 ft 


DO 7 1=74,165 
T 1 = 389.99 
PI = 464.86 
HI = 36500 
H = 500* (1-1 ) 

T(I) = T1 

P(I) = P1*EXP( -G0/ (R*T1 )* (H-H1 ) ) 
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RHO(I) = P(I)/(T1*R) 
7 CONTINUE 


Atmosphere from 82500 -> 100000 ft 


DO 10 1=166,201 
HI = 82500 
PI = 55.837 
T1 = 390.24 
A = . 001 646 
H = 500* (1-1 ) 

T(I) = (H-H1 )*A+T1 

P(I) = P1*EXP(-((G0/(R*T(I)))*(H-H1 ))) 
RHO(I) = P(I)/(T(I)*R) 

10 CONTINUE 


Dynanic pressures at various altitudes and mach nurbers 
Q{al t ,M) 


DO 15 ALTSUB= 1,201 
DO 20 MSUB=1 ,101 
M = 0.1* (MSUB-1 ) 

Q ( ALTSUB , MSUB ) = P(ALTSUB) * (GAM/2)*M*M 
20 CONTINUE 
15 CONTINUE 


WRITE ( 6 , * ) (’ Span (ft):’ ) 

READ(4, * ) SPAN 

WRITE ( 6 , * ) (’ Wing ref. area SE (ft‘2):’) 

READ(4, * ) SE 

WRITE ( 6 , * ) ( ’ Wetted wing area (ft*2) ') 

READ(4, * ) SAETW 

WRITE ( 6 , * ) (’ Wetted body area (ft *2):’) 

READ(4 , * ) S^ETB 

WRITE ( 6 , * ) (’ Max wing length (ft):’) 

READ( 4 , * ) 

WRITE ( 6 , * ) (' Body length (ft):') 

READ(4, * ) BL 

WRITE ( 6 , * ) (' Sweep angle (deg):’) 

READ(4 , * ) SAP 

WRITE ( 6 , * ) (' Max wing thickness (% chord):') 
READ(4 , * ) TC 

WRITE ( 6 , * ) (' Leading edge radius (% chord):’) 
READ(4,*) RLE 

WRITE ( 6 , * ) (’ Body cross-section area (ft A 2):’) 
READ(4 , * ) SB 

WRITE ( 6 , * ) (’ Body width (ft):’) 

READ<4, * ) BD 



WRITE ( 6 , * ) (' Starting Mach nurber (5 min):’) 

READ(5, * ) MS 

WRITE ( 6 , * ) (’ Ending mach nurber {10 max):’) 

READ(5, * ) ME 

WRITE ( 6 , * ) (’ Step size (.1 multiples):’) 

READ(5,*) MST 

WRITE(6,*) (’ Starting altitude (0ft min):’) 

READ (5 * ) ALS 

WRITE(6,*) (’ Ending altitude (100000ft max):’) 

READ(5, * ) ALE 

WRITE ( 6 , * ) (’ Step size (500ft multiples):’) 

READ(5, * ) ALST 

WRITE ( 6 , * ) (' Starting angle of attack:’) 

READ(5, * ) APS 

WRITE ( 6 , * ) (’ Ending angle of attack:’) 

READ(5, * ) APE 

WRITE ( 6 , * ) (' Step size:’) 

READ(5 f * ) APST 

WRITE ( 6 , * ) r Max drag:’) 

READ(5, * ) MAXD 

WRITE ( 6 , * ) (’ Max lift:’) 

READ(5, * ) MAXL 

WRITE < 6 , * ) (’ Min lift:') 

READ(5, * ) MINL 

DO 100 MACH=MS , ME , MST 

DO 110 ALT=ALS, ALE, ALST 

DO 120 ALPHAD=APS, APE, APST 

ALTSUB = ALT/ 500+1 
MSUB = MACH/. 1+1 
ALPHA = ALPHAD* ( PI / 1 80 ) 

MU = . 31 7E- 10* (T (ALTSUB ) **1 . 5* (734.7/ (T( ALTSUB )+21 6) ) ) 
CL1 = 2*<ALPHA/TC)*2/ (MACH* TC) 

CL2 » 2* (ALPHA/TC) * (2. 4/6*MACH*TC) * (4+( ALPHA/TC) **2 ) 

CL = (CL1 +CL2 ) *TC**2 

CD1 = 4/ (MACH-TC) * (4/3+ (ALPHA/TC)* *2) 

CD2 = 2. 4/3*MACH*TC* ( (ALPHA/TC) **4+8* (ALPHA/TC ) **2+16/5 ) 
CD = (CD1+CD2) *TC**3 

REB = (RHO( ALTSUB )*BL*MACH*SQRT(GAM*R*T( ALTSUB) ) ) /MU 
CF = .455/ (LOG10(REB) )**2.58 
CFCCFI = 1 / (1+. 144*MACH**2)** .65 
CDFB = CFCCFI *CF*SWETB/SE 

REW = (RHO( ALTSUB )*1AL*MACH*SQRT(GAM*R*T( ALTSUB) ) ) /MU 
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205 

200 

210 

215 


120 

110 

100 


CDFW = CFCCFI*CF*SWETW/ SE 
CD = COCDFW-CDFB 

LIFT = Q ( ALTSUB , MSUB ) *CL*SE 
DRAG = Q( ALTSUB, MSUB >*CD*SE 
LD = CL/CD 

IF( (LIFT.GT.MINL) .AND. ( LIFT . LT . MAXL ) .AND. 
MDRAG.LT.MAXD) ) TT-EN 

WRITE ( 1,205) MACH, ALT, ALPHAD 
FORMATdX, • MACH: ' ,F7.3, ’ ALT:',F12.4, 

*’ ALPHA: ’ ,F7.3> 

WRITE ( 1 , 200 ) LIFT, DRAG 

FORMATdX,’ LIFT=’,F12.5,’ DRAG=’ ,F12.5) 

WRITE(1,210) CL, CD 

FORMATdX,’ CL= ' , F7 . 4 , ’ CI>',F7.4) 

WRITE ( 1 ,215) Q( ALTSUB, MSUB), LD 
FORMATdX,’ Q=’ , F12.6, ’ L/D:’,F7.4) 

WRITE ( 2 , * ) CL 
WRITE ( 3 , * ) CD 
WRITE ( 7 , * ) MACH 
WRITE(8,*) ALT 
WRITE ( 9 , * ) LD 
WRITE (10,*) DRAG 
WRITEdl,*) LIFT 
WRITE ( 12, * ) ALPHAD 
ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 


STOP 

END 



Appendix D 

Inlet Analysis Program 


8 

1 

1 

26 

0 

380 
50000 
2.0 
4.0 
1 . 4 


Input 

nose angle 

num of cases 

number of outside ramps 

first inlet ramp 

number of inside ramps 

Required Mass flow rate 

altitude 

Mach No. 

alpha(angle of attack) 
gamma 


program inlet , . . . . _ . 

dimension pr ( 8 ) , rh ( 8 ) , te ( 8 ) , p0 ( 8 ) , rd{ 8 ) , areas { 8 ) , betas ( 8 ) 

dimension vs ( 8 ) , thetas ( 8 ) 
real mach,ml ,m2 ,mln,m2n 
integer stop 

common gam, theta, beta, ml , min, m 2 n , prat ,p0r, rrat , trat ,m2 

open ( l , f i 1 e= ' ' ) 
open( 2 , f i 1 e= ' ' ) 
pi = 4 . 0*atan (1.0) 
read( i,*)thetl 
read( 1 , * )numcases 
write( 2,450) 
do 900 k=l,numcases 
read( 1 , * )nout 
do 10 i =1 , nout 
read ( l, * ) thetas ( i+l) 

10 continue 

read( 1 , * )nin 
do 12 i=l,nin 

read( l ,*) thetas ( i+nout+1) 

12 continue 

thetas(nin+nout+2)=0.0 

read( l , * )wlr 

read( l , * ) al t 

read( i,*)mach 

read( l , * ) alpha 

read( 1 , * ) gamma 

gam = gamma 

thetas ( l ) = alpha+thetl 


20 


22 

23 


theta = thetas(l) 
ml = mach 
call oblique 
betas( l)=beta 
p0 ( 1 ) = p0r 
pr ( 1 ) = prat 
rh(l) = rrat 
te ( l ) = trat 
vs(l) = m2 
do 20 j=l,nout 
ml = m2 

theta = thetas ( j + l ) 
call oblique 
betas ( j + l ) * beta 
vs ( j+l ) = m2 
p0 ( j+l ) = p0( j ) *p0r 
pr ( j+l ) = pr ( j ) *prat 
rh( j+l ) = rh(j)*rrat 
te ( j+l ) = te( j ) *trat 
continue 
stop=0 


do 22 j=l,nin 
ml » m2 

theta = thetas( j+nout+1) 
call oblique 
if (m2. It. 1.0) then 
betas ( j+nout+1 ) = 0.0 
m2 = vs(j+nout) 


vs ( j+nout+1 ) 
p0( j+nout+1 ) 
pr ( j +nout + l ) 
rh( j+nout+1 ) 
te( j+nout+1 ) 
stop= j 
go to 23 
endi f 

betas ( j+nout+1 ) 
vs(j+nout+l) = 
p0(j+nout+l) = 
pr(j+nout+l) = 
rh(j+nout+l) = 
te( j+nout+l) = 
continue 
ml = m2 
call normal 
betas ( nout+nin+2 ) 
vs(nin+nout+2) = 
p0 ( nin+nout+2 ) = 
pr(nin+nout+2) = 
rh ( nin+nout+2 ) = 
te ( nin+nout+2 ) = 
call cond(alt, pres, temp, dens) 


vs ( j+nout ) 
p0( j+nout ) 
pr ( j+nout ) 
rh ( j+nout ) 
te( j+nout ) 


= beta 
m2 

p0( j+nout ) *p0r 
pr ( j+nout ) *prat 
rh( j+nout ) +rrat 
te ( j+nout ) *trat 


= 90. 
m2 

p0 ( j+nout ) *p0r 
pr( j+nout )*prat 
rh( j+nout )*rrat 
te ( j+nout ) *trat 
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pi = pr(nin+nout+2)*pres 

tl = te(nin+nout+2 ) *temp 

rl = rh(nin+nout+2)*dens 

ul = vs ( nin+nout+2 ) *sqrt ( gamma*1716*t 1 ) 

al = wlr*(pl/tl)*(518. 69/2116. 2)/rl/ul/32. 2 

areas ( nin+nout+2 ) = al 

sbmt =sin(pi*(betas( j+nout ) -thetas( j+nout) )/180. ) 
areas ( j+nout ) = al/sbmt 
do 30 i=l , j -2 

spb = sin(pi-pi*betas(nout+j-i+l)/180. ) 

sbmt=sin( pi * (betas ( nout+ j -i } -thetas ( j+nout -i ) )/180. ) 

areas ( nout+ j - i ) =spb*areas ( j+nout- i +1 ) /sbmt 

sbb = sin(pi*(betas( j+nout-i+l)-betas( j+nout-i ) 

$ +thetas( j+nout-i ) )/180) 

30 rd(nout+j-i )=sbb*areas( j+nout-i+l)/sbmt 

sb = sin(pi*betas(nout+2)/180. ) 

sbmt = sin(pi*(betas(nout+l)-thetas(nout+l) )/180. ) 
areas ( nout+1 ) =sb*areas ( nout+2 ) /sbmt 
spb = sin(pi-pi*(l. +betas ( nout+2 ) +betas (nout+1 ) - 
$ thetas(nout+l) )/180. ) 

rd ( nout+1 ) =spb*areas ( nout+2 ) /sbmt 
do 32 i=l,nout-l 

spb = sin(pi-pi*(betas(nout-i+2) )/180. ) 

sbmt=sin (pi* (betas (nout-i+1 ) -thetas ( nout - i +1 ) )/l80. ) 

areas ( nout- i +1 ) =spb*areas ( nout- i +2 ) /sbmt 

sbb = s in (pi *( betas ( nout- i+2) -betas ( nout -i+l)+ 

$ thetas(nout-i+l) )/l80. ) 

32 rd ( nout - i + 1 ) =sbb*areas ( nout -i+2 ) /sbmt 

sb = sin( be tas(2)*p i/180. ) 
areas(l) = areas(2)*sb 
ac = areas(l) 

print *, 'Capture area = ',ac,' nr = ' ,p0 ( nin+nout+2 ) 

wr i te ( 2 , 498 ) alt, mach, alpha, gamma 

write( 2, 500) ' M= ', (vs(j) ,j=l, nin+nout+2) 

wri te ( 2 , 500 ) 'P/Pinf = ',( pr ( j ), j=l , nin+nout+2 ) 

wr i te( 2 , 500 ) 'R/Rinf = ',( rh( j ), j=l, nin+nout+2 ) 

write( 2,500) 'T/Tinf = ',{ te( j ), j=l, nin+nout+2) 

write(2,500) 'Thetas = ',( thetas ( j ), j=l , nin+nout+2 ) 

wr i te ( 2 , 500 ) 'Betas = ', (betas ( j ), j=l , nin+nout+2 ) 

write( 2, 500) 'areas = ',( areas ( j ), j=l , nin+nout+2 ) 

wr i te ( 2 , 500 ) 'length = ',( rd( j ), j=l , nin+nout+2 ) 

wri te( 2 , 502 ) ac , p0 (nin+nout+2) 

900 continue 

450 format ( 20x, 'Mach 0-6 Inlet analysis') 

498 f ormat ( / , ' Al t i tude = ',f7.0,' ft. Mach No. - ', 

$ f 6 . 2 , ' Alpha = ' , f 6 . 2 , ' degrees Gamma = ',f5.3) 

500 format( lx,a9,8f 11.4) 

502 format ( 5x ,' Capture Area = ',f8.4,' Pr . ratio(nr) - ',fl0.6) 
stop 
end 
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subroutine oblique 

real thetbm,ml,mln,m2n,m2 

common gam, theta, beta, ml , min, m2n, prat ,p0r,rrat , t rat , m2 

pi = 4 . 0*atan (1.0) 

beta = thetbm(gam, theta, ml) 

min = ml*sin(beta*pi/180.0) 

rrat * (( gam+1 ) *mln**2 )/(( gam-1 ) *mln**2+2 ) 

prat = l + 2 *gam*(mln** 2 -i)/(gam+l) 

trat = prat/rrat 

m2n = sqrt ( (mln**2 + 2 . 0/ (gam-1 ) )/(2*gam*mln**2/( gam-l) -l)) 
m2 = m2n/sin( (beta-theta)*pi/180.0) 

rat={ 1 . 0+{ gam-1 . 0) *m2n**2/2 . 0 ) / ( 1 . 0+(gam-l . 0) *mln**2/2 . 0 ) 
p0r=prat* ( rat ) ** ( gam/ ( gam-1 . 0) ) 

print *, 'oblique =' ,ml ,mm, rrat ,prat , trat ,m2n ,m2 , p0r 

return 

end 


subroutine normal 
real ml ,mln,m2n,m2 

common gam, theta , bet a, ml , min ,m2n , prat , p0r , rrat, trat, m2 

m2 = sqrt ( ( 1 + ( gam- l ) *ml** 2 / 2 ) / ( gam*ml**2- ( gam- 1 ) /2 ) ) 

rrat = (gam+l)*ml**2/(2.0+(gam-l)*ml**2) 

prat = 1 + 2 . 0*gam* (ml**2-l)/ ( gam+1 ) 

trat = prat/rrat 

powg = gam/ ( gam- 1 . 0) 

rat=( 1 .+(gam-l. )*m2**2/2. )**powg/( l.+(gam-l . )*ml**2/2. )**powg 
p0r=prat* ( rat ) 

print *, 'normal =', ml, m2, rrat, prat, trat, p0r 

return 

end 


function thetbm( gam, th,ml ) 
real ml 

ffn(be)=(ml*sin(be) ) **2 - 1.0 

f f (be)=2*( f fn(be)/(ml**2*(gam+cos(2.0*be) )+2.0) )*cotan(be) 
pi - 0*atan (1.0) 
tth = tan(th*pi/l80.0) 
bl = 1.0 

1 tthl = ff(bl) 

if (abs(tthl-tth) . It. 0.000001) go to 999 
i f ( tthl . gt . 1 . 0 ) then 
bl = (pi/4.0+bl/2.0) 
go to 1 
endi f 

b2 = bl*0 . 8 

2 tth2 = ff(b2) 

i f ( abs (tth2-tth) . lt.0. 000001 ) then 
bl = b2 


go to 999 
endi f 

if (tth2.gt. l .0) then 
b2 = (bl+b2)/2.0 
go to 2 
endi f 

3 bnew = bl- (bl-b2 ) * ( tthl-tth) / ( tthl-tth2 ) 
if (bnew.gt. (pi/2.0) ) then 
bl = -1.0 
go to 999 
endi f 

tthnew = ff(bnew) 

i f ( abs ( tthnew- tth ) . 1 t . 0 . 000001 ) then 
bl = bnew 
go to 999 
endi f 
bl = b2 
tthl = tth2 
b2 = bnew 
tth2 = tthnew 
go to 3 

999 thetbm = bl*180/pi 

print * , ' thetbm = ' ,gam, th,ml, ' ->' ,bl*l90/pi 

return 

end 


subrout ine cond( al t ,pres , temp, dens ) 
if (alt. gt. 155349. 5)then 
temp = 508.79 

pres = 2. 4566*exp(-32. 2/1716/508. 79* ( al t - 156000 . 0) ) 
elseif (alt. gt. 82345.0) then 

temp = 389.99+0. 001631* (alt-82345.0) 

pres = 51.592 * ( temp/390 . 24 )**( -32 . 2/0 . 001631/ 17 16 ) 
elseif (alt .gt. 36 152.0) then 
temp = 389.99 

pres = 464. 86*exp(-32. 2/1716/389. 99* ( al t-36500 . 0) ) 
else 

temp = 518. 69-0. 00356*alt 

pres = 2116 . 2* ( temp/518 . 69 )**( 32 . 2/0 . 00356/1716 ) 
endi f 

dens = pres/1716/temp 

print *,'cond =', alt, pres, temp, dens 

return 

end 


